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ABSTRACT 

An  interactive  computer  code  for  the  preliminary  design 
of  solid  propellant  rocket  motors  ("SPRMD")  was  successfully 
developed  and  its  use  was  demonstrated  through  a  design 
example.  "SPRMD"  was  written  in  FORTRAN  for  use  on  an  IBM 
PC/AT.  It  combined  several  existing  codes  ("MICROPEP", 
"GRAINS",  "ROCKET",  etc.)  and  used  the  performance  loss 
estimation  methods  suggested  by  the  AGARD  Propulsion  and 
Energetics  Panel  for  aluminized  propellants. 
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I.   INTRODUCTION 

A.   BACKGROUND 

Solid  propellant  rocket  motor  design  is  a  very  complex 
process  which  often  requires  many  iterations  in  order  to 
accomplish  a  successful  design.  Usually,  the  first  step  in 
the  design  process  it  to  obtain  accurate  input  data  on  mis- 
sile system  restraints.  Examples  of  these  restraints  are 
ballistic  requirements  (which  can  include  thrust,  burning 
time,  total  impulse)  and  mission/vehicle  requirements  (which 
can  include  motor  length,  diameter,  allowed  hazard  and  plume 
properties,  altitude,  etc.).  These  parameters  are  often  re- 
ferred to  as  independent  parameters.  Through  the  use  of 
these  independent  parameters,  the  required  propellant  burn 
time  and  properties  can  be  roughly  estimated.  Then,  theo- 
retical motor  performance  and  losses  can  be  estimated.  The 
thrust-time  requirements  together  with  the  expected  struc- 
tural loads  permit  initial  selection  of  a  grain  design. 

The  next  step  of  solid  propellant  rocket  motor  design  is 
to  use  the  propellant  properties  and  grain  design  to  preuict 
the  pressure-  and  thrust-time  profiles  and  the  total  impulse 
of  the  motor.  Propellant  properties  and  grain  design  are 
then  iterated  until  these  independent  variable  requirements 
are  met,  together  with  acceptable  propellant  stress  levels. 
A  final   step  in  the   missile  design   process  is  to   use  the 
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thrust-time  profile  in  an  external  aerodynamics  code  to 
ensure  that  range  and  terminal  velocity  requirements  are 
met . 

During  the  past  decade,  there  has  been  significant  pro- 
gress in  the  automation  of  solid  rocket  motor  design.  It  is 
mainly  because  of  the  extensive  use  of  digital  computers, 
which  greatly  facilitate  the  iterative  type  computations 
mentioned  above  (Ref.  1).  "In  recent  years,  many  of  the 
design  processes  have  been  accomplished  using  separate  com- 
puter codes,  i.e.,  adiabatic  equilibrium  combustion  codes, 
codes  which  estimate  performance  losses,  grain  geometry 
optimization  codes,  etc.  These  codes  have,  in  some  cases, 
been  combined  into  much  larger  design  codes,  and  they  are 
often  difficult  and  expensive  to  run.  Recently,  PC  versions 
of  some  rocket  motor  design  codes  have  become  available. 
However,  they  must  still  be  used  individually,  and  were  not 
all  available  for  use  on  the  same  computer.  The  use  of 
these  codes  was  also  very  inconvenient  in  that  the  same 
parameters  often  had  different  variable  names  in  the  various 
codes.  The  user  had  to  transfer  or  convert  the  data  from 
one  code  to  the  next. 

In  addition,  a  recent  AGARD  publication  has  provided 
recomm.ended  emirical  equations  which  can  be  used  for  rapid 
and  accurate  estimation  of  performance  losses  for  aluminized 
propellants . 
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B.  OBJECTIVES 

The  first  objective  was  to  outline  a  suitable  "Solid 
Propellant  Rocket  Motor  Preliminary  Design  Procedure"  which, 
of  course,  is  not  a  unique  procedure.  The  second  objective 
was  to  create,  modify,  and/or  combine  the  PC  versions  of  the 
necessary  design  codes  into  one  interactive  code  which  could 
be  conveniently  used  for  preliminary  design  on  an  IBM  PC/AT 
microcomputer.  The  final  objective  was  to  demonstrate  the 
program  by  designing  a  tactical  solid  propellant  rocket 
motor  combustor. 

C.  APPROACH 

The  approach  to  achieve  the  objectives  was  first  to 
detail  an  acceptable  preliminary  design  process.  The  algo- 
rithm of  the  design  process  was  then  followed,  utilizing 
several  existing  codes  and  loss  estimation  methods  as 
appropriate,  to  obtain  an  integrated  design  code  for  use  on 
the  IBM  PC/AT.  The  programs  to  be  used  were  (1)  "MICROPEP" 
(a  microcomputer  version  of  a  propellant  evaluation  program 
from  China  Lake,  Naval  Weapons  Center),  which  calculates  the 
adiabatic,  equilibrium,  combustion  process,  and  theoretical 
motor  performance;  (2)  "GRAINS"  (a  radial  spoked  star  grain 
design  program  from  CSD) ;  (3)  "ROCKET  (a  Lockheed  internal 
ballistics  program) ,  which  is  a  one-dimensional  flow  code 
which  predicts  pressure  and  thrust  versus  time;  and  (4) 
"FLYIT"  (a  flight  simulation  program  from  China   Lake,  Naval 


12 


Weapons  Center)/  which  perforins  simplified  trajectory  calcu- 
lations. In  addition,  a  code  based  upon  AGARD-AR-230  was 
written  for  estimation  of  performance  losses  for  aluminized 
propellants.  Capabilities  of  the  final  program  were  then  to 
be  demonstrated  by  considering  the  design  of  a  tactical 
motor. 
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II .   SOLID  PROPELLANT  ROCKET  MOTOR  DESIGN  PROCESS 

A.  GENERAL 

The  method  of  solid  rocket  motor  design  adopted  in  this 

work  is  shown  in  Figure  1.   It  was  assumed  that  the  required 

information   on  fundamental   propellant   properties  (burning 
rates,   temperature    sensitivity,   etc.)    was    known.   In 

addition,   the   design   process   discussed   herein  did   not 

include   any   stress  analysis   for  the   propellant  grain  or 
case . 

B.  INPUTS 

The  inputs   that  are  required  to  satisfy   mission  objec- 
tives (the  independent  parameters)  were  as  follows: 

1 .  Ballistics  Performance  Requirements 

a.  Average  thrust  (F)  or  thrust-time  profiles, 

b.  Burn  time  (t^),  or 

c.  Total  impulse  (It^* 

These  parameters  are  interrelated  as  shown  in 

Figure  2.   Thus   It  =  J   ^^t  (2.1) 

o 

F  =  It/tb  (2.2) 

2 .  Mission/Vehicle  Constraints 

These  requirements  are  directly  related  to  the  motor 
design  and  include: 
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a.  Envelope 

Tne  allowable  envelope  (motor  length  L,  motor 
inside  diameter  D)  is  the  fundamental  constraint  on  the 
motor  geometry,  which  can  be  .  used  to  estimate  the  required 
burning  rate  (r)  and  total  volume  available  (V^)  for  the 
propellant  by  following  relations: 

r  =  (D/4)/tb  (2.3) 

(i.e.,  a  web  fraction  of  approximately  .5) 

Va  =  (  ttDV4)  L  (2.4) 

b.  Maximum  Expected  Operating  Pressure  (MEOP) 

The  MEOP  is  set  by  the  structural  limit  of  the 
motor  case.  The  maximum  operating  pressure  at  maximum  pro- 
pellant temperature  (Pmaxtp^  ^^  less  than  the  MEOP  by  the 
factor  of  reproducibility  tolerance  (see  Figure  2).  Pmaxtp 
can  be  initially  estimated  from  the  following  equation 

^maxtp"^'^^'-'^  ( 1-reproducibility  tolerance)         (2.5) 

c.  Temperature  Range  and  Altitude 

The  operating  propellant  temperature  range 
effects  the  selection  of  the  propellant  through  the  required 
burning  rate  temperature  sensitivity.  The  design  altitude, 
together  with  the  nominal  chamber  pressure  determine  the 
exhaust  nozzle  area  ratio. 
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d.   Allowed  Hazard  Rating  and  Allowed  Exhaust 
Signature 

The   hazard  rating  requirements  have  to   do  with 

whether   the   propellant   is   detonable    or   non-detonable . 

Modern  high   energy  propellants   are  often   detonable.    The 

allowed  exhaust   signature  affects  the   propellant  selection 

through  requirements  on  exhaust  smoke  and  radiation. 

C.   PROPELLANT  SELECTION 

From  the  ballistics  performance  and  mission/vehicle  con- 
straints, it  is  possible  to  make  a  preliminary  selection  of 
the  propellant  type  that  is  most  likely  to  give  the  required 
performance,  internal  ballistics,  flame  temperature  and 
mechanical  properties,  as  well  as  the  necessary  storage 
stability,  the  best  hazard  properties  and  exhaust  signature 
(Ref.  2).  A  particular  propellant  will  have  known  (or 
estimated)  burning  rate  (r),  temperature  sensitivity  (ttk)/ 
density  (pp)r   and  composition.    Then   the   nominal  chamber 

pressure   (Pc    )   can  be  calculated   by  using  the   equation 

nom 

( Ref .  3 ) : 


Pc     =  Pmaxtp  exP  (ttk  (Tnom  "  T^ax))  (2.6) 

nom 


D.   THEORETICAL  PERFORMANCE  COMPUTATION 

Following  the  propellant  selection,  theoretical  per- 
formance computations  can  be  made  for  the  specified  chamber 
and  ambient  pressures.    A  program  such  as  PEPCODE  (Ref.  4), 
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or  the  PC  version  "niCROPEP",  will  provide  data  such  as 
thrust  coefficient,  characteristic  exhaust  velocity,  shift- 
ing and  frozen  equilibrium  performance  results,  optimum  noz- 
zle expansion  ratio,  mole  fraction  of  condensed  species  in 
the  exhaust,  etc. 

E.  NOZZLE  GEOMETPxY  AND  MATERIAL  SELECTION 

The  information  on  theoretical  combustion  temperature, 
pressure,  gas  composition,  and  burn  time,  together  with 
geometry  constraints,  permits  the  nozzle  configuration  and 
material  to  be  selected. 

F.  PERFORMANCE  LOSS  ESTIMATION 

It  was  assumed  that  the  delivered  specific  impulse  could 
be  written: 

Icn  ~  Iqt-v    Hr    nr*  (2.7) 


sp  -  -^sp^^  ^Cf   ^C 


where   Ion.     is  the  theoretical  specific  impulse  calculated 
^Pfh 

by  "MICROPEP";  and  nc  ^^^  ^C*  ^^®  ^^^  thrust  coefficient 
efficiency  and  characteristic  velocity  efficiency  (or  com- 
bustion efficiency),  respectively. 

nc    ^^s   computed   by   summing   the   effects   from   the 
following  losses 

1.  Divergence  losses  (eoiv^' 

2.  Two-phase  flow  losses  (  e^j-p )  , 

3.  Boundary  layer  loss  (ebl^' 
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4.  Kinetics  loss  (ekin^' 

5.  Submergence  loss  (esUB^' 

6.  Erosion  loss  ( eerOS ^  • 

Thus,  TIC   =  1~(  ^DIV"^^TP"^^BL'*'^KIN"*"^SUB"^^EROS)/100        (2.8) 

The  combusion  efficiency  depends  primarily  on  the 
residence  time  (tj^)/  and  the  latter  can  be  expressed  as 

tj.  =  (Vol/m)  (P^  M/12  R  T^)  (2.9) 

where 

Vol  i   ^he  instantaneous  volume  of  the  combusor  (in^), 

p^   is  the  chamber  pressure  (Ibf/in^), 

T^   is  the  chamber  temperature  ( °R) , 

M    is  the  molecular  weight, 

-  R    is  the  universal  gas  constant  (ft  Ibf /Ibj^Qj^^  "^R)  ,and 

-  m    is  the  mass  discharge  rate  (Ibm/sec). 

For  example,  Figure  3  (Ref.  2)  is  a  plot  of  specific  impulse 
versus  average  residence  time  for  aluminized  propellants. 
It  shows  that  the  residence  time  must  be  somewhere  between 
10  to  15  msec  to  obtain  a  combustion  efficiency  between  90% 
and  95%.  Increasing  the  free  volume  of  the  combustion 
chamber  will  improve  the  residence  time.  In  this  work  the 
initial  volume  was  used  in  the  residence  time  calculation, 
i.e.,  the  worst  case. 

G.   THROAT  AREA  AND  EXIT  AREA 

Accurately  sizing  the  throat  and  exit  areas  is  one  of 
the  major  goals  of  the  solid  propellant  rocket  motor  design. 
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FIGURE  3, 
EFFECT  OF  RESIDENCE  TIME  ON  ^sp  EFFICIENCY (Ref  2) 
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These  areas  are  dependent  on  several  other  variables, 
such  as  thrust  coefficient  efficiency,  combustion 
efficiency,  specific  impulse,  and  etc..  In  the  sizing 
process,  it  is  convenient  to  start  by  using  the  theoretical 
values.  The  effects  of  losses  can  then  be  handled  itera- 
tively  until  the  solution  converges  to  the  final  areas.  The 
rest  of  this  section  will  present  the  basic  algorithm  for 
the  throat  and  exit  area  sizing  process. 

The  basic  definition  of  specific  impulse  can  be  ex- 
pressed as  (Ref.  3): 

IsD    =  f^/^  (2.10) 

^^th 

where 

F  =  Cf    PcAth  (2.11) 

th 


Then 


^fth  "^Cf  Pc 


The  theoretical   throat  area   can  be   estimated  by  using 

Cf     for   shifting  eguilibrium   flow  (from  'MICROPEP"  -  the 
th 

theoretical  performance  computation).  The  losses  which  are 
calculated  to  ootain  nc  depend  upon  the  output  from 
"MICROPEP"  and  the  throat  area.  Thus,  iteration  is  required 
to  obtain  both  the  losses  and  throat  area. 

In  summary,   the   calculation   can  be   conducted  in   the 
following  steps: 
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1.  Calculate  A^.^  using  equation   (2.11). 

2.  Calculate  nc  ^^  discussed  above. 

3.  Calculate  A|-h  using  equation  (2.12). 

4.  Repeat   steps  2.  and  3.  until  the  values  for  throat 

area  and  nr   converge, 
^f 

In  the   current  program   the  iteration   was  terminated   when 
there  was  less  than  a  1%  change  in  A^j^  P®^  iteration. 
The  exit  area  can  then  be  calculated  from 
Ae  =  A^h  e  (2.13) 

where  e  is  the  exit  to  throat  area  ratio.  e  is  calculated 
in  "MICROPEP",  assuming  that  the  nozzle  exit  pressure  is 
equal  to  the  local  ambient  pressure  (one-dimensional,  ideal 
flow  is  assumed). 

H.   GRAIN  DESIGN 

The  primary  purpose  of  grain  design  is  to  shape  the  pro- 
pellant  grain  such  that  the  rocket  motor  can  deliver  the 
required  performance  (thrust-time  or  pressure-time  profile). 

In  order  to  be  able  to  select  the  general  configuration 
of  the  grain  geometry,  several  parameters  must  first  be 
determined:  the  web  fraction  (Wf),  volumetric  loading  frac- 
tion (V^) ,  and  the  length-to-diameter  ratio  (L/D)  (Ref.  5). 
L/D  is  known  from  the  input  data.  Web  fraction  is  the  frac- 
tion of  the  motor  radius  that  is  filled  with  propellant. 
Thus, 


23 


The  propellant  burning  rate  (r)  is  generally  expressed  in 

the  form: 

n 
r  =  apc  (2.15) 

where  a  and  n  are  empirical  constants.  They  are  assumed  to 
be  known  for  the  selected  propellant.  The  volumetric  load- 
ing fraction  is  defined  as  the  ratio  of  the  volume  occupied 
by  the  propellant  (V^)  to  the  total  volume  available  in  the 
motor  for  propellant  (V^).   Thus, 

Vj;  =  -T-2-  (2.16) 


where 


Va  =  (^)  L  (2.17) 


The  required  propellant  volume  is  found  as  shown  below 

It  =  Ftb  =  Ispi^tb  =  Isp  r^  tb  =  IspWp  (2.18) 


Thus,  Wp  =  It/Isp  (2.19) 

This  yields  the  required  propellant  weight  (Wp)  and 


Vd  =   -^  (2.20) 


Pp 

where  pp  is  the  density  or  specific  weight  of  the  selected 
propellant.  With  Wf ,  VJl,  and  L/D  known,  an  initial  selec- 
tion of  grain  geometry  can  be  made  (for  example,  see  Ref. 
5).  For  a  selected  grain  configuration,  the  web  versus 
burning   surface    area   (which   is   the   data    needed   for 
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thrust   and   pressure   versus   time   calculations)    can   be 
computed . 

I.   THRUST  AND  PRESSURE  VERSUS  TIME  CALCULATIONS 

The  thrust  and  pressure  versus  time  calculations  are 
made  using  conservation  of  mass  for  the  combustor (Ref .  3). 
The  mass  addition  from  the  burning  propellant  (Wg)  is 
egual  to  the  sum  of  the  mass  accumulated  in  the  combustor 
and  the  mass  exhausted  from  the  nozzle  (va^)  as  shown  in 
Figure  4.   Thus, 

irig  =  dM/dt  +111^  (2.21) 

mrj  is  the  mass  of  gas  generated  from  the  burning  propellant. 


g  =  iPp^h^   =  ZPpAbaPc  (2.22) 


The  summation  is  over  all  propellant  grains  that  are  burning 

at  the  same  time. 

dM    .    ^  -,   ■    .    , 

-rr        IS  the  mass  accumulating  in  the  combustor. 

dM     d_,   „  ,.      dVol  ^  ,,  T  dpg 

where       .  (2.23) 

dVol    r-      ,„^    ^Pg  ^   1   dPc 
__  =  IrAi,        and   -^    -.  -— 


1 


then 


dH4        Pc      dVol         vol      dP 


c 


^^         RT.       ^^  RT  '^t 


+   - —      —rr  (2.24) 
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FIGURE    4. 


MASS    CONSERVATION    IN    A    SOLID    PROPELLANT    MOTOR 
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nii^  is  the  mass  exhausted  from  the  nozzle.   For  chocked  flow. 


iTin  - 


Pc  ^th  gc 


(2.25) 


Combining  equations  2.21  through  2.25  results  in 

V     A   =,0"   -       (       ^C^   dVol,  ,V01^^PC     ,    r^C^th^Ci  (J       ^f-s 

1  RTc  c 

Solving  for  dP^/^t 

^=  Ol^^ca  PpA^aP?  -  ^4i^)  -  Pe(^)]    (2.27, 

Equation  2.27  can   be  numerically   integrated  to  give  P(^(t). 
The  thrust-time  behavior  can  then  be  obtained  from 

F  (t)  =  nc   Cf (t)Pc(t)Ath(t)  (2.28) 

J.   TERMINAL  VELOCITY  AND  RANGE  CALCULATIONS 

The  last  step  of  solid  propellant  rocket  motor  design  is 
to  incorporate  the  designed  motor  into  the  desired  airframe 
and  make  a  flight  simulation.  The  purpose  of  this  flight 
simulation  is  to  verify  that  the  motor  can  produce  the 
desired  range  and  terminal  velocity. 
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III.   MICROPEP 

"MICROPEP"  is  a  microcomputer  version  of  the  propellant 
Evaluation  Program  (PEP),  published  in  Naval  Weapons  Center 
Report  NVJC-TP-6037 ,  "Theoretical  Computations  of  Equilibrium 
Compositions,  Thermodynamic  Properties,  and  Performance 
Characteristics  of  Propellant  Systems,"  by  D.  R.  Cruise, 
April  1979  (Ref.  4).  The  code  (MICROPEP)  can  handle  12  dif- 
ferent input  ingredients  and  200  combustion  species.  It 
performs  the  calculations  for  the  thermodynamic  properties 
and  performance  of  propellant  systems. 
A.   BASIC  ASSUMPTIONS 

The  basic  assumptions  of  "MICROPEP"  are  as  follows: 

1.  Equilibrium  adiabatic  combustion  at  the  nozzle 
entrance , 

2.  One-dimensional,  isentropic  expansion  in  the  nozzle  in 
which  the  flow  can  be  either  a  shifting  equilibrium 
flow  (assuming  chemical  equilibrium  is  maintained 
throughout  the  expansion  process)  or  frozen  composi- 
tion flow  (assuming  the  chemical  composition  of  the 
the  flow  remains  the  same  throughout  the  expansion 
process ) . 

3.  Gases  were  considered  to  be  perfect. 

4.  Nozzle  exit  pressure  is  assumed  to  be  equal  to  the 
ambient  pressure. 
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B.  INPUTS 

Inputs  required  by  "MICROPEP'  are: 

1.  Composition  and  mass  of  each  ingredient.  Normally  the 
weights  of  ingredients  are  chosen  to  add  up  to  100  q. 
NOTE:   Ingredients   may  be   input  by   the  user  or  they 

may  be  selected  from  an  existing  data  file.   In 

the  latter   case  ingredient   serial  numbers  are 

used. 

For  assistance,  use  README,   which  should  be  on 

the  PEPCODED.DAF  Data  File. 

2.  Combustion  pressure. 

3.  Nozzle  exit  pressure.  This  defines  the  limits  for  the 
expansion  process. 

C.  OUTPUTS 

The  program  calculates: 

1.  Combustion  product  composition,  properties  ,  and 
temperature  at  the  nozzle  entrance. 

2.  Nozzle  exit  product  composition,  properties,  and 
temperature  for  two  cases: 

a.  Shifting  equilibrium  flow 

b.  Frozen  composition  flow. 

3.  Required  throat  area  and  exit  area  with  no  losses,  for 
the  same  two  cases. 

4.  Theoretical  specific  impulse,  thrust  coefficient,  and 
characteristic  exhaust  velocity  for  both  shifting 
equilibrium  and  frozen  composition  flow. 
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IV.   LOSS  MECHANISMS 

For  estimation  of  the  performance  losses,  a  set  of 
empirical  formulas  suggested  by  the  Propulsion  and 
Energetics  Panel  of  AGARD  in  AGARD-AR-230  (Ref.  6)  were 
used.  Similar  methods  can  also  be  found  in  NASA  Report  SP- 
8039  (Ref.  7)  "Solid  Rocket  Motor  Performance  Analysis  and 
Prediction."  It  is  assumed,  as  in  equation  2.7,  that  the 
combustion  efficiency  (nc*)  ^'^^  thrust  coefficient  effici- 
ency (  ric  )  c^i^  be  treated  ^  independently.  The  empirical 
formulas  are  valid  only  for  aluminized  propellants. 

A.   THRUST  COEFFICIENT  EFFICIENCY 

The  thrust  coefficient  efficiency  (  ric  )  ^^^  ^^  expressed 
as  in  equation  2.8 

^C   =  l-(  SDIV+eTP+£BL'*"^KIN'''eSUB"*"eEROS)/100*0 

Empirical  equations  for  each  of  the  losses  (e's)  will  be 
presented  and  briefly  discussed  below.  More  details  can  be 
found  in  Reference  6. 

1 .   Divergence  Loss  (^DIV) 

This  is   the  loss   due  to  the   radial  velocity   com- 
ponent of  the  gas  at  the  nozzle  exit  plane 


eoiv  =  50  [1-Cos  i^-^r^)]  .  (4.1) 
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a  is  the  nozzle  half  angle  and  6ex  is  the  exit  angle  for  a 
contoured  nozzle.  For  a  contoured  nozzle,  a  is  measured  to 
a  line  drawn  from  the  exit  and  tangent  to  the  wall  at  the 
throat . 

2,  Kinetics  Loss  (^KIN) 

This  is  the  loss  (or  reduction  in  performance  from 
the  value  of  shifting  equilibrium  flow)  associated  with  not 
attaining  chemical  equilibrium  throughout  the  nozzle  expan- 
sion process.  It  is  estimated  to  be  1/3  of  the  fractional 
difference  between  the  specific  impulse  for  shifting  equi- 
librium (Isp  )  ^"^d  frozen  composition  (Igp  ).  Thus, 

£KIN  =  33.3  [1  -  -^^  ]  (4.2) 

-^sps 

Both  ■'■sps  and  ^spf  are  obtained  from  'MICROPEP." 

3 .  Boundary  Layer  Loss  ( ^BL) 

The  boundary  layer  loss  (ebl^  ^^  expressed  as 


p  0.8  p  0.8 

EBL  =  Ci   ^  Q  ^  [1  +  2  exp(-C2   ^  ^  ^  t^)]  [1  +  0 .  016 ( e-9 ) ] 

Dth  *  D^h 

(4.3) 


This  loss  accounts  for  both  reduced  flow  areas  and  transient 
heat  loss.  In  equation  4.3,  P^  is  in  psi,  D^.^^  is  in  inches. 
The  C]^  and  C2  coefficients  are 
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for  ordinary  nozzles         C-^   =  0.00365 

C2  =  0.000937 
for  steel  nozzles  Ci  =  0.00506 

C2  =  0.0 
4 .   Two-Phase  Flow  Loss  ( ^TP ) 

This  loss  is  due  primarily  to  the  velocity  lag  of 
the  condensed  species  of  the  combustion  products  as  the 
mixture  passes  through  the  exhaust  nozzle. 


where 


MfC4     D^^S 


a.  D^^  is  the  diameter  of  throat  in  inches 

b.  Dp  is  the  mean  AI2O3  particle  diameter  in 

microns  and  is  calculated  from 
Dp  =  3.39  Dth°-^^^2 

c.  P(2    ^^    chamber  pressure  in  psia 

d.  C's  are  dependent  on  the  diameter  of  the  throat 
C4   =  0.5 

Dth   <  1:  C3  =  9  C5  =  1  Ce  =  1 

1  <  Dth  <  2:       C3  =  9  C5  =  1  Cg  =  0.08 

Dth>2  and  Dp<4:    C3  =  13.4  C5  =  0.8  Cg  =  0.8 
Dth>2  and 

4<Dp<8:  C3  =  10.2  C5  =  0.8  Cg  =  0.4 

Dth>2  and  Dp>8:    C3  =  7.58  C5  =  0.8  Cg  =  0.33 

e.  Mf  is   the  mole  fraction  of  consdensed  phase  in 
moles/100  gm  of  reactants. 
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5 .   Submergence  Loss  (^SUB) 
This  loss  is  expressed  as 


PcE^O.8    qO.^ 
AtTT 


esuB  =  0-0684  (— ^)     — ^^-^  (4.5) 


Dth 
It  is  due  to  the  nozzle  configuration  for  a  submerged  nozzle 

where  S  is  the  submergence  length  in  inches  (the  length  that 

the  nozzle  is  imbedded  in  the  combustion  chamber). 

6 .   Nozzle  Erosion  Loss  ( ^EROS ) 

This  loss   results  from   nozzle  erosion   during   the 

burn  and  is  expressed  as 

^EROS  =  tl  -  r|§^l  ^  100  (4.6) 

Ic-,    is  the  theoretical  I,,^  for   the  initial   nozzle  expan- 

sion  ratio   and  ^sp-^    is   the   theoretical  Igp  for   the   mean 

expansion  ratio.    ^^Pm      ^^    determined   by  interpolation   of 

lcr>  versus   the    nozzle   expansion    ratio   data    from 

^Pth 

"MICROPEP." 

B.   COMBUSTION  EFFICIENCY 

NASA  Report  SP-8064  (Ref.  2)  states  that  the  combustion 
efficienty  (nc*)  is  determined  by  "the  completeness  of  metal 
combustion  within  the  motor  and  by  the  degree  to  which  com- 
bustion products  reach  chemical  equilibrium  among 
themselves."  These  effects  depend  primarily  on  the 
residence  time  (t^-)  as  expressed  in  equation  2.9.   An  nc*  of 
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90-95%   can  be  attained  if  residence  time  is  10-15  msec  (see 
Figure  3 ) . 

C .   PROGRAM 

The  subroutines  which  solve  for  the  thrust  and  combus- 
tion efficiencies  are  listed  in  Appendix  A,  Subroutine  PFE 
and  Subroutine  LOSCF.  Basically,  they  followed  the  method 
discussed  in  Chapter  II.   The  algorithm  is  shown  below 

1.  Initialize  Cf  (coefficient  of  thrust)  as  ^f^h 
(theoretical  coefficient  of  thru Jt  from  "MICROPEP" 
for  shifting  equilibrium) . 

2.  Calculate  A^h  ^^i^  ^th* 

Dth  =  2(     p  J 
^fth  c^ 

3.  Calculate  nr-   using  the  relation 

^f 

i^C^  =  l-(  eDlv'*"SBL+£TP+^KIN+^SUB"*"SEROs)/100 

4.  Calculate  Cf  using 

Cf  =  Cf    nc 
th    f 

5.  Correct    D^j-^   using 

Dth   =    2    (-i^ ) 

CfP^TT 

6.  Repeat  Steps  3  to  5  until  A^h  converges. 

The  reason  that  A^h  ^^^  used  for  the  iteration  variable  was 
that  He   (through  e3L,  Etp  '  ^'^^    ^SUB^  depends  upon   A^j^* 
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V.   GRAINS 

"GRAINS"  is  PC  based  code  for  the  design  of  radial  sided 
star  grains.  The  original  version  was  written  by  G.  J. 
V7oten  and  R.  McCormick  of  CSD,  United  Technologies,  November 
1984.  It  was  then  modified  by  LT  G.  Liston  of  AFWAL/Port, 
December  1984.  The  general  shape  of  this  grain  and  the 
geometric  variables  are  shown  in  Figure  5. 

"GRAINS"   can  be   used   in   either   of  two   modes:    the 
"Design  Mode"  or  the  "Burnback  Model." 
A.   "DESIGN  MODE" 

In  the  "Design  mode"  the  program  solves  for  geometries 
that  satisfy  input  reguirements  on  the  number  of  star  points 
reguired  port  area,  perimeter  factor,  web,  and  etc., 
together  with  the  allowed  tolerances  on  these  parameters. 

The  geometric  constraints  for  the  grain  are  determined 
from  the  propellant  characteristic  and  motor  constraints. 
For  example: 


Aport  =  (|)dM(1  -  (j-^)] 


D,  t]3  and  F  are  motor  inputs  reguirements  and  I^p  and  p^  are 
determined  once  the  propellant  has  been  selected. 
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FIGURE  5. 


NOMENCLATURE  FOR  "GRAINS"  DESIGN  MODE 
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B.  "BURNBACK  MODE" 

The   "burnback   mode"   uses   the  length   of  the   grain, 

initial  port   area,  perimeter,  and  the  results  from  "Design 

Mode"   to  calculate   the  perimeter,   perimeter  factor,  port 

area,  and   burning   surface  area   as  functions  of   the   web 
burned . 

C.  COMMENTS  ON  PROGRAM  OPERATION 

Care  must  be  used  for  the  inputs  and  input  tolerances, 
otherwise  no  acceptable  solutions  will  be  obtained.  It  is 
best  to  start  with  large  allowable  tolerances  on  ApQj-^,  etc. 
(typically  20%).  After  initial  solutions  are  found  in 
"Design"  for  the  number  of  star  points,  range  of  web 
thicknesses,  etc.,  the  range  of  the  input  variables  can  be 
reduced  along  with  the  allowable  tolerances. 
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VI.   ROCKET 

"ROCKET"  is  also  a  PC-based  computer  code.  It  was 
originally  written  by  The  Lockheed  Company  for  a  mainframe 
computer.  Later,  it  was  transferred  for  use  on  the  IBM-PC 
by  J.  P.  Francis  of  the  Naval  Weapons  Center  (August  1985). 
The  function  of  "ROCKET"  is  to  calculate  the  pressure-time 
and  thrust-time  profiles  for  the  motor,  including  theoreti- 
cal calculations  for  both  ignition  and  tailoff  transients. 
The  description  of  the  motor  can  include  igniters,  multiple 
propellant  grains,  and  motor  insulation.  In  this  thesis 
work  the  program  was  restricted  to  a  motor  with  a  single 
grain  configuration. 
A.   BASIC  ASSUMPTIONS 

The  assumptions  made  in  "ROCKET"  are: 

1.  The  combustion  products   are  ideal  gases. 

2.  Propellant   burning   rate   follows    the   expression 

n 
r  =  aPc  . 

3.  There   are  no   effects  of   mass  addition   or  erosive 
burning  in  the  combustion  chamber. 

4.  C*  varies  with  P^  in  the  form  of 


C*  =  C*ef(^^) 


^C   >iX 

Pref 
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5.  The  temperature  sensitivity  (ttj^)  is  used  to  correct 
for  the  effects  of  propellant  temperature  on  the 
propellant  burning  rate  and  chamber  pressure. 

r        =  r.vo^)  e0.01a^(Tp-70°) 
P 

where 

Qp  =   (  1-n)  TTj^ 

B.  ANALYTICAL  BACKGROUND 

The  program  is  based  upon  the  method  discussed  in 
Chapter  II,  i.e.,  it  numerically  integrates  eguation  (2.27) 
to  obtain  P^(t).  Then  it  calculates  C^  with  a  fixed  ratio 
of  specific  heats  (y).  Thrust  is  then  calculated  using 
equation 

F  =  Cf (t)Pc(t)Ath(t) XCq 
where   X   is  nozzle   divergence  loss   and   Cq   is  the  nozzle 
discharge  coefficient.    X  and  Cq  ^^^    only  two  parts  of  n^  . 
In  order   to  use  ROCKET  in  the  combined  program,   X  was   set 
to  unity  and  C^  was  equated  to  riQ  • 

C.  INPUTS 

Inputs  to  "ROCKET"  include: 
1.   General  motor  parameters: 

a.  Throat  area  (A^h^ 

b.  Nozzle  exit  area  (Aq) 

c.  Total  motor  volume  (V^) 

d.  Ambient  temperature  and  pressure 
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e.  Nozzle  throat  plug  closure  blowout  pressure 
(nozzle  is  pluqged  until  this  pressure  is 
exceeded,  typically  35  psia) 

f.  Initial  pressure  in  the  motor  (Pzero'  typically 
15  psia) 

g.  Nozzle  divergence  loss  X  (set  equal  to  1  in  this 
work ) 

h.   Ratio   of  specific   heats   for   chamber   gases  (y) 
(this    was   obtained    from    "MICROPEP"   chamber 
conditions ) 
i.   Design  pressure  (P^) 

j.   Throat  radial  erosion  rate   (erosr). 
2.   Motor  element  descriptions: 

a.  Propellant  burning  rate  (r) 

b.  Burning  rate  exponent  (n) 

c.  Propellant  temperature  sensitive  (ttk) 

d.  Burning  rate  reference  pressure 

e.  Propellant  C* 

f.  Pressure  correction  exponent  for  C* 

g.  Propellant  density  (pp) 

h.   A  table  of  web  versus  burning  surface  area. 

D.   OUTPUTS 

The  outputs  include: 

1 .  A  pressure  versus  time  table 

2.  A  thrust  versus  time  table 
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3.   The  free   volume,   propellant   consumed   (weight),  and 
throat  area  versus  time  table. 
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VII.   FLYIT 

"FLYIT"  is  a  three-degree-of-f reedom  (range,  altitude 
and  angle  of  attack)  trajectory  simulation  program  origi- 
nally developed  by  G.  Burgner  of  the  Propulsion  Analysis 
Branch,  Naval  Weapons  Center  (NWC).  A  User's  guide  was 
later  published  by  Y.G.  Coenen  (Ref.  8).  The  program  can 
simulate  vertical-plane  trajectories  which  include  take-off, 
climbs,  cruises,  landing,  and  rocket  boost.  It  is  a  tool  to 
synthesize  and  analyze  the  trajectories  of  air-to-surface, 
and  surface-to-surface  missiles. 

The  basic   structure  of   "FLYIT"  consists  of   five  func- 
tional sub-groups,  they  are: 
Atmospherics 
Trajectory  controls 
-   Aerodynamics 
Rocket  boosters 
Sustanier  propulsion. 
Users  are  allowed  to  select  a  different   option   in  each 
functional  group  to  simulate  different  flight   profiles  with 
different  propulsion  systems.  Through  the  use  of   "FLYIT"  it 
is  possible  to   estimate  the  terminal  velocity  and   range  of 
certain  propulsion  systems. 

"FLYIT"  was  not  incorporated  into  the  main  body  of  the 
present  program  due  to  time  limitations. 
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VIII.   COMBINED  PROGRAM 

The  combined  program  ("SPRMD")  is  listed  in  Appendix 
and  basically  follows  the  design  process  discussed  in 
Chapter  II.  The  basic  structure  of  this  program  is  shown  on 
Figure  6. 

This  program  combines  "MICROPEP"  (to  calculate  the 
theoretical  thermodynamic  performance  of  the  propellant), 
"PFE,"  and  "LOSCF"  (to  calculate  ncf/  ^"^  to  estimate  nc*f 
and  size  the  throat  and  exit  areas),  "GRAINS"  (to  provide 
web  burned  versus  burning  surface  area  for  radially  sided 
star  grains),  and  "ROCKET"  (to  compute  the  pressure-time 
and  thrust-time  profiles).  Inputs  to  the  program  include 
"INPUTO"  (the  initial  mission  requirement  inputs),  "INPUTl" 
(inputs  for  "MICROPEP"),  "INPUT2"  (inputs  for  "GRAINS"),  and 
"RCKTIN"  (inputs  for  "ROCKET").  All  of  these  programs  com- 
prise an  interactive  program  for  solid  propellant  rocket 
motor  preliminary  design. 

A.   ALGORITHM 

A  simplified   algorithm  of   this  program   can  be  briefly 
listed  as  follows: 
1.   "INPUTO" 

Input   the    basic    ballistic    and    mission/vehicle 
requirements.      Calculate    the    required    average 
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INPUT   VARIABLES 


"MICROPEP"     AT  2     P's 


250  psi,     Pc 


THEORETICAL  PERFORMANCE 


AGARD  LOSS  CALCULATIONS 


-> 


V 


ESTIMATED  ACTUAL  PERFORMANCE 


•GRAINS" 


± 


Ab  VS.  WEB  BURNED 


"ROCKET' 


THRUST  VS.  TIME 
PRESSURE  VS.  TIME 


\/ 


"FLYIT3" 


_-^._. 


TERMINAL  VELOCITY 
RANGE 


FIGURE  6. 


FLOW  CHART  FOR  THE  COMBINED  PROGRAM 
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requirements.  Calculate  the  required  averaqe 
chamber  pressure  for  a  specific  propellant  and 
propellant  temperature. 

2.  "INPUTl" 

Input  propellant  composition,  chamber  pressure,  and 
ambient  pressure.   Create  a  file  for  'MICROPEP." 

3.  "MICROPEP" 

Perform  the  theoretical  performance  calculations  for 
the  selected  propellant. 

4.  "PFE" 

Performs  the  following  functions: 

a.  Estimates  the  combustion  efficiency 

b.  Calculates  the  delivered  specific  impulse 

c.  Sizes  the  throat  area,  exist  area,  calculates 
perimeter  factor,  perimeter,  required  web, 
required  burning  rate,  and  etc. 

5.  "INPUT2" 

Input  the  geometrical  description  of  the  star  grain 
and  propellant  characteristics.  Create  a  file  to 
accommodate  the  processing  of  "GRAINS."  Other  grain 
design  programs  (such  as  SPP)  can  be  readily 
incorporated. 
6.  "GRAINS" 

Designs  the  star  grain  and  calculates  web  versus 
burning  surface  area. 
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7.  "RCKTIN" 

Uses  the  results  from  the  previous   programs  to  create 
a  file   for  "ROCKET." 

8.  "ROCKET" 

Calculates  the  pressure-time  and  thrust-time  profiles. 

B.   BRANCHING 

The  algorithm  discussed  above  also  provides  for  the 
capability  of  branching.  For  example,  if  the  desired  P^ 
versus  time  is  not  attained,  the  user  can  either  re-do  the 
grain  design  and/or  select  other  propellant  properties. 
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IX.   EXAMPLE  OF  DESIGN  PROGRAM 

The  initial  requirements  which  are  input  into  the  pro- 
gram are  shown  in  Table  1,  together  with  a  selected  value 
for  the  temperature  sensitivity  of  the  propellant  (irj^).  A 
reproducibility  tolerance  of  7%  was  chosen. 

The  nominal  operating  pressure  at  Tp  =  70°F  was  calcu- 
lated to  be  1030  psia  and  the  required  propellant  burning 
rate  for  a  web  fraction  of  0.5  was  0.615  in/sec. 

A  composite  propellant  consisting  of  78%  AP ,  12%  HTPB, 
and  10%  Al  was  selected.  "MICROPEP"  was  then  run  for  pres- 
sures of  250  psia  and  the  nominal  P^  value  of  1030  psi. 
Outputs  from  "MICROPEP"  are  presented  in  Table  2. 
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TABLE  1 
INPUT  PARAMETERS 


AVERAGE  THRUST,  LBF 14250 , 0 

BURN  TIME,  SEC 6.50 

MINIMUM  THRUST  REQUIRED,  LBF 8000.000000 

MOTOR  LENGTH,  INCH 40.00 

MOTOR  DIAMETER,  INCH'.  . 16.00 

MAXIMUM  EXPECTED  OPERATING  PRESSURE...  13  00.0 

TPMAX,  F. 150.0 

TPNOM,  F 70.0 

TPMIN,  F -65.0 

TEMPERATURE  SENSITIVITY,  %/F. 0.200 

NOMINAL  CHAMBER  PRESSURE ,  PSI 1030.2 

REQUIRED  BURNING  RATE,  IN/SEC 0.  615 

ALTITUDE,  FT 10000. 

AMBIENT  PRESSURE,  PSI 10.13 

AMBIENT  TEMPERATURE,  R 483  . 1 


'i 


48 


TABLE  2 


MICROPEP  OUTPUT  FOR 


P^  =  250  AND  1030  PSIA: 


P^  =  10.0  PSIA.  78%  AP,  12%  HTPB,  10%  Al 


DESIGN  1 


Case 


1  of 


DH 


DENS  COMPOSITION 


AMMONIUM  PERCHLORATE  (AP) 

HTPB  (SINCLAIR) 

ALUMINUM  (PURE  CRYSTALINE) 


-602   0.07040     ICL   4H    IN 
13   0.03320   103H   73C    10 
0   0.09760     lAL 


40 


INGREDIENT  WEIGHTS  (IN  ORDER)  AND  TOTAL  WEIGHT 
78.0000     12.0000     10.0000    100.0000 


(LAST  ITEM  IN  LIST) 


INGREDIENT  VOLUME  RATIOS 

AMMONIUM  PERCHLORATE  (AP)       7  0.487% 

ALUMimm  (PURE  CRYSTALINE)       6.518% 


HTPB  (SINCLAIR) 


THE  PROPELLENT  DENSITY  IS  0.063  62  LB/CU-IN  OR   1.7610  GM/CC 

NUMBER  OF  GRAM  ATOMS  OF  EACH  ELEMENT  PRESENT  IN  INGREDIENTS 

3.895570  H      0.878965  C      0.663847  N      2.667428  0 
0.370645  AL     0.663847  CL 
CP  -       43.62668610 

CHAMBER  RESULTS  FOLLOW: 

T(K)   T(F)   P(ATM)    P(PSI)  ENTHALPY   ENTROPY    CP/CV      GAS 
3248.  5388.    17.01    250.00    -46.80    240.56   1.1990   3.644 

DAMPED  AND  UNDAMPED  SPEED  OF  SOUND=  3211.743  AND  3563.987  FT/SEC 

HEAT  CONTENT  (298  REF)   1253.397  CAL/GR       2257.631  BTU/LBM 


RT/V 
4.667 


22.995% 


MOL  WT 
27.44 


0.98321  H20 
0.32835  N2 
0.07546  HO 
5.93E-03  02 
1.32E-04  A1C13 
2.31E-05  CHO 
7.13E-06  NH3 
2.03E-06  CNH 


0.71991 
0.17947 
0.07476 
4.51E-03 
.28E-04 
,73E-05 
,55E-06 
,91E-06 


1. 
1. 
5. 
1. 


CO 

A1203* 

CI 

AlOCl 

C12 

NH 

NHO" 

A1H02 


CP 


41.92294693 


EXHAUST  RESULTS  FOLLOW: 

T(K)   T(F)   P(ATM)    P(PSI)  ENTHALPY 
2313.  3704.     0.69     10.13   -108.91 


0.58380  H2 
0.15900  C02 
0.00880  O 
3.65E-03  A1C12 
1.19E-04  AlO 
1.61E-05  N 
5.24E-06  A120 
1.09E-06  CH20 


3 

ENTROPY 
240.56 


CP/CV 
1.1995 


0.57339  HCl 
0.11256  H 
0.00709  NO 
3.20E-03  AlCl 
6.37E-05  Al 
1.21E-05  H02 
4.85E-06  AlH 


GAS 
3.510 


RT/V 
0.196 


DAMPED  AND  UNDAMPED  SPEED  OF  S0UND=  2650.051  AND  2952.060  FT/SEC 
HEAT  CONTENT  (298  REF)    851.166  CAL/GR        1533.130  BTU/LBM 


MOL  WT 
28.49 


1.01958  H20 
0.33181  N2 
0.01421  CI 
1.33E-04  A1C12 
4.03E-05  AlCl 


0.67108  CO 
0.20785  C02 
0.00497  HO 

1.04E-04  O 

1.25E-05  A1C13 


0.64916  HCl 
0.18284  A1203& 
0.00233  A1203* 
l.OlE-04  AlOCl 
9.06E-06  C12 


0.59218  H2 
0.01785  H 
0.00021  NO 
7.51E-05  02 
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TABLE  2  CONTINUED) 


MICROPEP  OUTPUT  FOR  Pc  =  250 


P^  =  10.0 


PSIA.  78%  AP,  12% 


AND  1030  PSIA: 
HTPB,  10%  Al 


PERFORMANCE:   FROZEN  ON  FIRST  LINE,  SHIFTING  ON 


ISP    GAMMA   T* 
226.7   1.1835  2975. 
232.5   1.1497  3045. 


P*     C* 
9.65  5060.1 
9.77  5082.8 


CF     OPT  EX 
1.442    4.31 
1.472    4.72 


SECOND  LINE 

D-ISP    A*M.  EX  T 

399.2  0.62924  1976. 

409.4  0.63205  2313. 


BOOST  VELOCITIES  FOR  PROPELLANT  DENSITY  OF  0.06362   (S.G.  OF    1.761) 

5./23451.  10./18585.  15./15857.  25./12612.  30./11520.  55./  8215. 
60./  7788.  69./  7128.  71./  6998.  88./  6064.  100./  5548.  150./  4113. 
175./3647.   '00V3277.   300./2336.  1000./  780.  3000./  269.  5000./  163. 


EXP. 

RATIO 

1. 
2. 
3. 
4. 

5. 

6. 

7. 

8. 

9. 
10. 
11. 
12. 
13. 
14. 
15. 
16. 
17. 
18. 
19. 
20. 
21. 
22. 
23. 
24. 
25. 


EXIT 
PRESS 
ATM 
9.769 
3.256 
1.215 
0.833 
0.644 
0.523 
0.439 
0.378 
331 
294 
264 
240 
219 
202 
137 
0.174 
0.163 
153 
144 
136 
129 
123 
117 
0.112 
0.107 


0. 
0. 
0. 
0. 
0. 
0. 
0. 


,9 
,0 
,4 
,3 
,0 
.5 
,3 


EXIT 

PRESS 
SI 

989.6 

329. 

123. 
84. 
65. 
53. 
44. 
38. 
33.5 
29.8 
26.8 
24.3 
22.2 
20.5 
19.0 
17.7 
16.5 
15.5 
14.6 
13.8 
13.1 
12.4 
11.8 

■  11.3 
10.8 


EXIT 

TEMP 

K 

3044. 

2678. 

2386. 

2314. 

2313. 

2311. 

2310. 

2310. 

2309. 

2308. 

2307. 

2307. 

2306. 

2306. 

2305. 

2305. 

2305. 

2304. 

2304. 

2304. 

2303. 

2303. 

2303. 

2302. 

2302. 


OPTIMUM 
IMPULSE 
SEC 
104.9 
175, 
214. 
226. 
234. 
240.7 
245.8 
250.1 
253.8 
257.1 
260.0 
262.6 
265.0 
267.2 
269.2 
271.1 
272.8 
274.4 
276.0 
277.4 
278.8 
280,1 
281.3 
282.5 
283.6 


OPTIMUM 

IMPULSE 

SI 

1029. 

1716. 

2103. 

2223. 

2300. 

2361. 

2410. 

2453. 

2489. 

2521. 

2550. 

2575. 

2599. 

2620. 

2640. 

2658. 

2675. 

2691. 

2706. 

2720. 

2734. 

2746. 

2758. 

2770. 

2781. 


VACUUM 

IMPULS 

SEC 

195.7 

235.5 

248.3 

257.6 

264.5 

269.9 

274.4 

278 

281 

284 

287 

289 

291 

293 

295 

297 

298 

300 

301 

302 

304 

305 

306 

307 

308 


VACUUM 

IMPULS 

SI 

1919. 

2309. 

2435. 

2526. 

2593. 

2647. 

2690. 

2728. 

2760. 

2789. 

2814. 

2837. 

2858. 

2878. 

2896. 

2912. 

2928. 

2942. 

2956. 

2969. 

2981. 

2992. 

3003. 

3014. 

3024. 


SEA  LV 
IMPULS 
SEC 
186.4 
216.9 
220.4 
220.5 
218.0 
214 
209 
203 
197 
191 
184.8 
177.9 
170 
163 
155 
148 
140 
132.8 
124.9 
116.9 
108.9 
100.8 
92.6 
84.4 
76.1 


SEA  LV 
IMPULS 
SI 
1846. 
2149. 
2184. 
2184. 
2160. 
2121. 
2074. 
2019. 
1960. 
1897. 
1831. 
1762. 
1691. 
1619. 
1545. 
1469. 
1393. 
1316. 
1237. 
1158. 
1079. 

998. 

917. 

836. 

754. 
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TABLE  2  (CONTINUED) 

MICROPEP  OUTPUT  FOR  P^  =  250  AND  1030  PSIA: 
p   =  10.0  PSIA.  78%  AP.  12%  HTPB.  10%  Al 


DESIGN  1    Case   2  of 


DH    DENS  COMPOSITION 


AMMONIUM  PERCHLORATE  (AP)         -602   0.07040     ICL   4H    IN    40 
HTPB  (SINCLAIR)  13   0.03320   103H   73C    10 

ALUMINUM  (PURE  CRYSTALINE)  0   0.09760     lAL 

INGREDIENT  WEIGHTS  (IN  ORDER)  AND  TOTAL  WEIGHT         (LAST  ITEM  IN  LIST) 

78.0000     12.0000     10.0000  \   100.0000 

INGREDIENT  VOLUME  RATIOS      •     '^ 

AMMONIUM  PERCHLORATE  (AP)     ,,70.487%    HTPB  (SINCLAIR)        .         22.995% 

ALUMINUM  (PURE  CRYSTALINE)     "  6.518% 

THE  PROPELLENT  DENSITY  IS  0.06362  LB/CU-IN  OR   1.7610  GM/CC 

NUMBER  OF  GRAM  ATOMS  OF  EACH  ELEMENT  PRESENT  IN  INGREDIENTS 

3.895570  H      0,878965  C      0.663847  N      2.667428  O 
0.370645  AL     0.663847  CL 
CP  =»       43.75085831 

CHAMBER  RESULTS  FOLLOW: 

T(K)  ■T(F)   P(ATM)    P(PSI)  ENTHALPY   ENTROPY    CP/CV      GAS    RT/V   MOL  WT 
3375.  5616.    70.08   1030.24    -46.80    230.37   1.1952   3.597    19.485    27.80 

DAMPED  AND  UNDAMPED  SPEED  OF  SOUND=-  3247.693  AND  3  603.495  FT/SEC 

HEAT  CONTENT  (298  REF)   1309.442  CAL/GR       2358.580  BTU/LBM 

1.01084  H20  0.71716  CO  0.59375  HCl  0.57367  H2 

0.32895  N2  0.17941  A1203*      0.16171  C02  0.07526  H 

0.05727  HO  0.05233  CI  0.00586  NO  0.00506  A1C12 

4.49E-03  O  3.70E-03  AlOCl  3.15E-03  02  2.63E-03  AlCl 

3.14E-04  A1C13  1.83E-04  C12  7.24E-05  AlO  5.77E-05  CHO 

3.63E-05  Al  2.75E-05  NH  2.71E-05  NH3  1.55E-05  N 

1.33E-05  H02  9.39E-06  NHO  8.06E-06  CNH  5.08E-06  AlH 

4.45E-06  CH20  3.50E-06  A120  2.16E-06  A1H02 
CP  =       40.18283081 

EXHAUST  RESULTS  FOLLOW: 
t 
T(K)   T(F)   P(ATM)    P(PSI)  ENTHALPY   ENTROPY    CP/CV      GAS    RT/V   MOL  WT 
1852.  2874.     0.69     10.13   -130.14    230.37   1.2087   3.491     0.197    28.64 

DAMPED  AND  UNDAMPED  SPEED  OF  SOUND=  2373.392  AND  2644.581  FT/SEC 

HEAT  CONTENT  (298  REF)    659.260  CAL/GR       1187.465  BTU/LBM 

0.98547  H20         0.66303  HCl         0.63198  CO  0.63022  H2 

0.33192  N2  0.24696  C02         0.18531  A1203&  0.00097  H 

7.87E-04  CI  l.lOE-04  HO  2.78E-06  A1C12  2.25E-06  NO 

2.21E-06  A1C13  1.52E-06  HH3 

PERFORMANCE:   FROZEN  ON  FIRST  LINE,  SHIFTING  ON        SECOND  LINE 
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TABLE  2  (CONTINUED) 


MICROPEP  OUTPUT  FOR  Pc  =  250 
P^  =  10.0  PSIA.  78%  AP,  12% 


AND  1030  PSIA: 
HTPB,  10%  Al 


ISP    GAMMA  T* 
262.5   1.1907  3081. 
269.3   1.1617  3139. 


P*     C* 
39.69  5084.9 
40.09  5126.9 


CF    OPT  EX 
1.661   12.28 
1.690   13.25 


D-ISP    A*M.    EX  T 
462.2  0.15344   1610. 
474.2  0.15471   1852. 


BOOST  VELOCITIES  FOR  PROPELIANT  DENSITY  OF  0.06362   (S.G.  OF    1.7^1) 


5./27165. 

10./21528.   15 

./18368. 

25./14609.   30./13344.   55, 

/  9516 

60./ 

9021. 

69./  8257.   71 

./  8106. 

88./  7024.  100./  6426.  150. 

/  4764 

175./4224. 

200./3796.   300./2705. 

1000./ 

904.  3000./  312.  5000./  188 

EXP. 

EXIT 

EXIT 

EXIT 

OPTIMUM 

OPTIMUM 

VACUUM 

VACUUM 

SEA  LV 

SEA  LV 

RATIO 

PRESS 

PRESS 

TEMP 

IMPULSE 

IMPULSE 

IMPULS 

IMPULS 

IMPULS 

IMPULS 

ATM 

SI 

K 

SEC 

SI 

SEC 

SI 

SEC 

SI 

1. 

40.059 

4057.9 

3139. 

106.6 

1045. 

197.7 

1939. 

195.4 

1936. 

2. 

13.353 

1352.6 

2722. 

176.8 

1734. 

237.6 

2330. 

233.0 

2308. 

3. 

4.904 

496.8 

2391. 

216.7 

2125. 

250.1 

2453. 

243.3 

2410. 

4. 

3.289 

333.2 

2270. 

229.5 

2250. 

259.4 

2544. 

250.3 

2479. 

5. 

2.435 

246.7 

2183. 

238.3 

2337. 

266.0 

2609. 

254.7 

2523. 

6. 

1.913 

193.8 

2115. 

245.0 

2402. 

271.1 

2558. 

257.4 

2550. 

7. 

1.563 

158.3 

2060. 

250.3 

2454. 

275,1 

2698. 

259.2 

2568. 

8. 

1.314 

133.1 

2014. 

254.6 

2497. 

278.5 

2731. 

260.3 

2579. 

9. 

1.129 

114.3 

1975. 

258.3 

2532. 

281.4 

2759. 

260.9 

2584. 

10. 

0.986 

99,9 

1940. 

261.4 

2563, 

283.8 

2783. 

261.1 

2586. 

11. 

0.873 

88.4 

1910. 

264.2 

2590. 

286.0 

2805. 

261.0 

2585. 

12. 

0.781 

79.1 

1882. 

266.6 

2615. 

287.9 

2824. 

260.7 

2582. 

13. 

0.706 

71.5 

1858. 

268.8 

2636. 

289.7 

2841. 

260.1 

2577. 

14. 

0.643 

65.1 

1835. 

270,8 

2655. 

291.3 

2856. 

259.4 

2570. 

15. 

0.589 

59.7 

1814. 

272.6 

2673. 

292.7 

2870. 

258.6 

2562. 

16. 

0.543 

55.0 

1795. 

274.3 

2690. 

294.0 

2883. 

257.6 

2552. 

17. 

0.503 

51.0 

1778. 

275.8 

2705. 

295.2 

2895. 

256.6 

2542, 

13. 

0.468 

47.4 

1761. 

277.2 

2718. 

296.4 

2906. 

255.5 

2531. 

19. 

0.438 

44.3 

1746. 

278.5 

2731. 

297.4 

2917. 

254.2 

2519. 

20. 

0.411 

41.6 

1731. 

279.8 

2743. 

298.4 

2927. 

253.0 

2506. 

21. 

0.386 

39.1 

1717. 

280.9 

2755. 

299.4 

2936. 

251.6 

2493. 

22. 

0.365 

36.9 

1704. 

282,0 

2766. 

300.3 

2944. 

250.2 

2479. 

23. 

0.345 

34.9 

1692. 

283.0 

2776. 

301.1 

•  2952. 

248.8 

2465. 

24. 

0.327 

•33.1 

1681. 

284.0 

2785. 

301,9 

2960. 

247.3 

2450. 

25. 

0.311 

31.5 

1670. 

284.9 

2794. 

302.6 

2968. 

245.8 

2435. 
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The  output  from  "MICROPEP"  indicated  the  following 

C*    =  5126.9  ft/sec 
th 

I-^   =  269.3  Ibf  sec/lbm 
^Pth 

e  =  13.25 

Moles  of  condensed  AI2O3/IOO  gm  =  0.1794 
Pp  =  .0636  Ibm/in^ 
The  AGARD-AR-230  empirical  eguations  were  then  used  to 
estimate  the  value  for  nc  •  '^c*  ^^^  chosen  to  be  0.93.  For 
^  the  calculated  residence  time  of  7.2  msec.  Figure  3,  indi- 
cates that  a  better  value  may  have  been  approximately  0.87. 
As  discussed  above,  the  losses  are  iterated  with  the  throat 
area.  The  resulting  losses  and  throat  area  are  presented  in 
Table  3. 

Once  the   losses  have   been  estimated   Igp  can  be  deter- 
mined from 

Isp  =  ^C^^C*  -"-sp  j^ 

The  following  parameters  were  then  calculated: 
m  =  F/Igp 

Vp  =  It/^spPp 
Vil  =  Vp/V3 


Dpe  =  (/T^V7)D 

Ap  =  (j^Dpe  =  I  d'  (1-V,) 


J    =  Ath/Ap 
These  parameters  are  shown  in  Table  4. 
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TABLE  3 
ESTIMATED  PERFORMANCE  LOSSES 


PERCENT  DIVERGENCE  LOSS 1 . 704 

PERCENT  TWO  PHASE  FLOW  LOSS 3.202 

PERCENT  BOUNDARY  LAYER  LOSS 1.828 

PERCENT  KINETIC  LOSS 0.849 

PERCENT  SUBMERGENCE  LOSS 0  .  000 

PERCENT  EROSION  LOSS 0.007 

THRUST  COEFFICIENT  EFFICIENCY 0.S24 

Etacf=        0.924   Cfc=        1.562   Ath=  8.856 

PERCENT  DIVERGENCE  LOSS 1.704 

PERCENT  TWO  PHASE  FLOW  LOSS 3  .  199 

PERCENT  BOUNDARY  LAYER  LOSS 1.783 

PERCENT  KINETIC  LOSS 0.849 

PERCENT  SUBMERGENCE  LOSS 0.000 

PERCENT  EROSION  LOSS 0  .  007 

THRUST  COEFFICIENT  EFFICIENCY 0.925 


54 


TABLE  4 
CALCULATED  MOTOR  DATA 


THRUST  COEFFICIENT .  1.  563 

ETACF 0.925 

SPECIFIC  IMPULSE,   LBFSEC/LBM 231.587 

MASS  FLOW  RATE,     LBM/SEC 61.532 

PROPELLANT  VOLUME,   CUIN 6286.8 

NOZZLE  THROAT  AREA,  SQIN 8.851 

NOZZLE  EXIT  AREA,    SQIN 117.290 

EQUIVALENT  PORT  DIAMETER,  IN 7  .,476 

PORT  AREA,  SQ  IN  43.893 

PERIMETER  FACTOR 2  .  489 

REQUIRED  WEB,  IN 2.500 

REQUIRED  WEB  FRACTION 0  .  312 

VOLUMETRIC  LOADING 0  .  782 

LENGTH  TO  DIAMETER  RATIO 2  .  500 

RESIDENCE  TIME,  MSEC 7.256 

THROAT  TO  PORT  AREA  RATIO  (J  FACTOR)...  0.202 

REQUIRED  BURNING  RATE,    IN/SEC 0.385 
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At  this  point  the  following  expressions  remain  to  be 
satisfied: 

Perimeter  Factor  =  PFAC  =  Pwetted^'^^pe 

'PFAC  =Ai3/TTDpgL 

^b  =  ^wetted  ^ 

m  =  Pp  A^  i:p 

^  cnom 

r  =  web/^5 

In  these  three  expressions  the  unspecified  variables  are 
PFAC,  Aj^,  r,    and  web. 

For  a  selected  web  or  r  the  other  parameters  can  be 
determined.  For  example,  in  the  current  problem  Table  E 
summarizes  the  required  values  for  selected  web 
thicknesses . 


TABLE 

5 

REQUIRED  PARAMETERS 

FOR  SELECTED 

WEBS 

web( in) 

r  (in/sec) 

Ak  (in2) 

PFAC 

3.5 

.538 

1796 

1.91 

3.0 

.461 

2095 

2.22 

2.5 

.385 

2514 

2.67 

Once  the  required  burning  rate  (r)  is  determined,  appro- 
priate   values   of  a    and    n    are    required   for   the 

n 
expression  r  =  ^P^^  . 

The  grain  design  problem   is   to   find   a  geometry  which 

will  satisfy  Table  5,  together  with   the  required  conditions 

on  L,  D,  PQ-t  and  F-t  (and  stress). 
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In  the  present  program/  only  radially  spoked  star  grains 
can  be  analyzed.  It  may  not  be  possible  to  satisfy  the 
motor  requirements  with  this  type  of  grain.  The  process  is 
essentially  one  of  trial-and-error . 

Three  attempts  at  meeting  the  design  requirements  are 
presented.   One  6-point  and  two  8-point  designs  were  tried. 

A  typical  output  from  the  "Design  flode"  of  "GRAINS"  is 
shown  in  Table  6. 

Once  the  "GRAINS  Design  Mode"  has  been  run,  "GRAINS 
Burnback"  is  run  to  obtain  A5  versus  web. 

The  output  from  "Burnback"  for  Case  II  is  shown  in  Table 
7. 

"ROCKET"  can  then  be  run  to  obtain  the  pressure-time  and 
thrust-time  profiles.  Input  data  for  Case  II  are  shown  in 
Table  8. 

A  small  section  of  the  output  from  "ROCKET"  is  shown  in 
Table  9  for  Case  II. 

A  summary  of  the  results  for  the  three  cases  is  pre- 
sented in  Table  10  and  pressure-time  traces  are  shown  in 
Figure  7. 

None  of  the  cases  attempted  met  the  desired  design 
values  with  sufficient  accuracy.  Further  iterations  may  re- 
sult in  a  better  solution.  However,  the  star  design  has 
considerable  sliver  with  the  resulting  long  tailoff.  This 
section  was  presented  in  order  to  demonstrate  the  process 
required  to  obtain  a  final  grain  design. 
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TABLE  6 
OUTPUT  FROM  "DESIGN  MODE"  OF  "GRAINS" 


Design  i.'iDUtS 

!1axiaUB  NuRber  of  SDokes  =  8     brain  Outer  Radius  =   S. 00000 

Desired  Port  Area  =  43.8'3293    Desireo  Peri»eter  Factor  =   2.48915 

RI!  =    3.50000  DEL2  =  -0.25000    FIN2  =  2.00000 

SRIl  =   0.25000  DELI  =  0.10000    FINl  =  0.25000 

SRI2  =   0.25000  DEL4  =  O.lOOOd    FlH4  =  0  25000 

FINWEB  =  4.00000  UEBSTP  =  0. 25000    KSFOKE  =  0.50000 

CIRTOL  =  0.10000  PTTOL  =  0.10000    -ERTQL  =  0. .10000 

Tolerarice  Lints 


For:  Area:  fliniaua  =   39.50364  ilaxiaua  =  ^3. 23222 

Periieter:  faniaui  =   52.61331  f.axiaui  =  64.30515 

Circle   :  fliniaua  =  324.00000  riaxiEua  =  336.00000 

lieo  Tnickness  =   2.5000 


'   j.^.i;"   ---'M     Ttr-             "^o*      01       •>")  1.-  'inc  c"*'  C?*^  Cl'r"  C'J'I  0-''!  SU'T    r,»T  • 

>  rtrvKi      -zr^in         rr-w         p»l<i         Ki           f.i  R*  r.Kj  ij^i  ,   bf.i  o.-Fm  ^rtil  rriii  rnu       "rti'* 

ERROR  ->  ttl  greater  tnan  PTTGL,  therefore  eaaal 

i   Di.b:^  ww.w.:w  /.lu-.   j.jUI'  j. 'jv   ...^jV  j.jv'J  cjUU  v...j»  ■.■.j.j'.'  j...  ^.s^-j  ^u.si-  .../^J  j..'s^ 

,8  ^:..3^-^  :t.:t?  2.3{4-  :.:!:   :.5:    :.:i,  :.5c:  i.o:-  >.::.  :.:!  :.i:-  4. 09o  :£.;.-  :.^::  -.:3:  iii 

9  42.:::  52.737  2.::;  :.w   ■.:5o  5.:5:  :.5i.';  ^.00;  :.:5c  ..:!:  :.::};  -.i::  :,.3ii  .,"::  ^.:::  11 


1--.U-  -.-■  ii^  cfsc^rr  tnan  -  iw_,  trereior-  eoui. 

s  ii.,ii/   a..i.ii  i.a.j  i../-v  j.'jv'J  j.ijj  -.juv  o.OOU  v...jv  v.^jU  j.i/-  t./oi,  ^o.o.t  „,ii.i    ....as 

ERROR    -;■    Ql  creaier  tnar;  -'TQ;.,  trierefore  eciial 

^  W...JC:  aj.Ov'..     j.:7w     ..j'JO     .:.rjl'     j.^jl-  j.jUv     o.U'Jv     O.^jU     0.<.Dv     j.J7;     j.j.io  ^c.i^i'T     i.,ti.j     i.,i.j^ 

ERROR  -.'  2i  Greater  ;rkan  PTTOL,  rnerefore  eaual 

ERROR  ->  S2  greater  tian  PERTOL,  therefore  ecual  -           .  '      ■ 

5  27.651  55.344  3.532  2.250  2.500  5.250  5.500  S.OOC  5.250  C.250  3.440  5. 733  26.814  2.723  1.24E 

ERROR  ->  rhi4:[  O.O:  less  tnan  0.0 

SpoKes:  6  Phil:  0.1  Phi2:  0.5  Phi3:  0.0 
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TABLE  6  (CONTINUED) 
OUTPUT  FROM  "DESIGN  MODE"  OF  "GRAINS" 


iieb  Thicicness  =  2.7500 


K  AFOR"  ?ERI«   'FAC   SRI   R2  R3    R4    RR5   .SRI   SRI   SPKT   PHI!  .PHI2   FHI3  ?HI4 

ERROR  ->  S2  crtd;?-  than  FERTGL,  tnc'efcre  ecual. 

TJifSJ     tSi>3.          ^■ViD         wlwVV/'          Mai'uV  waUUv  w  t  ^UV          wa'v'VV          Va^JV          V  •  i^w'v          wli.0/'        Os^.-J     Jv  i  /  w  3          ^>Uw::          'JtSi* 

ERROR  -)  CI  create'  tnan  PERTOw,  tne^'efcre  *Quai 

-1  '  •■,  •1"';  c<^  "'(■c   "  ir  ■   -■  -.r.'-,   "^  'C,',  ^    •'^•■"""  -  '-'^   3  ''sf '■   '■  ""^-"^     ""^i"'   ■<■'-■   i  (-j.- -  "  ■  ■"<'■•  r    ■-'-   •■»  ""^ 

C   T'^/t^wki  Jv's  r  Vw    *.K..w'T    w«4^tiV    wvJVv  w.vV^  Mi«.Uv    StVVv    Jt^JV    .  ■  ^w^    Waw^^    ^iV^C   Mv'irC*    ^bCOD    V  ■  ^0  / 

ERROR  ->  >^ni4:;  0.0^  Usi  tnan  0.0 

Ssokes:  3  =n:l:  0.1  Phi2:  0.5  Phi3:  0.1 
ERROR  -;  rhi4:i  O.t]   ies5  tnan  O.C 

ScoKis:  B  Phil:  O.l  ?'i:2:  0.5  'h:2:  0.1 
ERROR  ->  ?r!i4::  O.OI  Uss  than  0.0 

isoKes:  B  Phil:  0.1  Phi2:  0.5  ?hi3:  0.1 
ERROR  -;  Phi-:1  -0.1]  Uss  than  0.0 

ERROR  ->  Phi4::  -O.i:  i^ss  tnar  CO 

Sookes:  3  Phil:  0.1  ?hi2:  0.5  ?n:3:  O.i 
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TABLE    7 
GRAINS     "BURNBACK"    OUTPUT    FOR    CASE    II 


*axii:-B  Nuaoer  or'  Spcses  =    0  jram  Oute''  Sacius  =       3.00000 

Desired  Port  Area  =      41. 12200         Desire:  Periae;er  Factor  =       2.5530C 


RIl  = 

3.00000 

DEL2  = 

0.00000 

riN2  = 

0.00000 

SRIl  = 

0.25000 

DELI  = 

0.00000 

riNi  = 

C. 25000 

S-I2  = 

0.25000 

DEL4  = 

0.00000 

r  Iri*  = 

C. 25000 

flKtita  - 

0.00000 

yEBSTP  = 

0.00000 

:<SPQ.<E  = 

0.00000 

tiRit'i.  = 

0.00000 

P-TCL  = 

0.00000 

PER'QL  = 

o.ocooo 

Burnoacic  Incuts  ^ 

^rh\r\  Lengt'!  =      40.00000         FHI2  =      2£.S1400         PHI4  =       :.rO3O0 

STEPl  =       5.         STEP2  =       5.         STEP3  =       Z.      '    STEP*  =       5.         STEPS 

NuiDer  of  Spokes  =8 

mat  Thickness  =       2.5C0O0 

Ueb       Periaeter  Permeter  Burn  Port 

Burned  Factor  Area  Area 

(IN)  (IN)  (3Q  IN)  (5u  ::, 


&urnback  of  S^I  (Spoke  Corner  Radius) 

0.000       58.737  2.553  2349.472  42.121 

0.05C       53.051  2.481  2352.033  45.065 

0.100       53.365  2.416  2374.604  48.027 

0.15C       5S.679  2.357  2387.170  51.003 

C.200       59.*993  2.303  2399.737  53.394 

0.250       60.308  2.253  2412.303  57.002 

0.250       60.308  1.253  2412.303  57.002 

61.389  2. 054  2455.576  71.054 

£2.150  1.904  2435.997  84.747 

62.589  1.784  2503.564  37.973 

62.707  1.682  2508.279  110.643 

61.504  1.591  25C:.141  122.639 
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TABLE    7     (CONTINUED) 
GRAINS    "BURNBACK"    OUTPUT    FOR    CASE    II 


1.323  c:.50i  :.57S  2500.141  124. 8St 

1.335  £0.310  1.511  2422.411  128.313 

1.4&S  53.117  1.447  2364.680  132.817 

1.540  57.424  :.2Bi  2236.350  136.533 

1.613  55.730  1.328  2229.220  140.247 

^•DOD  J-r.Jw;  >.../.  ..d.tar  ItC/i/' 

cjrnDick  of  refsBinoer  to  nil  -  C 

1     '7^  c     t;;;  ■      ()Q  I've;    ttn  'C'     '02 

"^  T:!!  ^;    q?A  >    ■■lyc  if'''   toe  '77   "i^ 

2.500  52.^4:  1.0S5  2037. S74  1=6. 12i 


Lsro  ^iti  ih;ci;n*55 


Percen:  of  fuel  rrsainin?  =      3. 40   , 


Burnnact!  of  reaaining  fuel 

2.832  2C.320         v. 403  812.733  136.073 

2.343  15.430         0.303  £17.134  133.051 

3.066  10.413         C.20S  416.753  :33.637 

3.200  5.27B         0.105  211.107  200. 6S3 

3.343  COCO         0.000         C.OOC  2C  1.055 
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TABLE    8 
INPUT    DATA    FOR    "ROCKET" ,    CASE    II 


lotcr  PerfcrBanc*  -rograh  l'i*-r'^  version  l.^; 

f«DTGR  HAS  1    e^AINS 

DESCRIPTION  CF  6RAIN  1 

BCRI^i  RATE  ^^'^^^^ 

BURN  RAT£  EXP.  ^.4000 

PI  SUB  K  ^-20^^ 

BURN  REr.  PRESS.  1030.2 

p  cTiC  4763.0 

C  STAR  EX?.  ^-^^^i 

0     1*  C  ' 

isNI'IQN  TIHE  O.OOOp 

DElTA  IcK.  tiwe  ^-^^^^ 

-  Noraal  Incut  15  finisned 

Seneral  ConfiQuration  ?ARAi*.ETEF.3  are  ... 


Tnroat  AREA 

5.3510 

Exit  AREA 

117.2903 

Esaansion  RATK! 

1  o      -VCI  T 

Iota;  flotor  VCLL'ME 

Ovti. J 

TeBOt-rarare 

/O.O 

AilBIEST 

►TeiSUre 

10.1 

Closure 

CUllOui 

35.0 

Pzero 

15. C 

Laabda 

1.0000 

C.324£ 

':a£&a 

1.1332 

Throat  1 

design  r'eSSiir 

*  v  w  V  .  fc 

THROAT  EROS.  DELAY  Tl     2000.000 
RESSURE  vs.  THROAT  ERSSIGN  RATE 


600. 

0.000000 

1400. 

0.000000 

^ 

UEB 

BURh  AREA 

0.000 

2343.472 

0.250 

2411.952 

I.IOB 

2507.744 

1.540 

2296.212 

2.011 

2060.360 

2.500 

2097.107 

2.B32 

812.473 

3.342 

O.OOO 

TOTAL  PROPELLANT  WEISHT= 


411.99b  I6N1TER  WE16HT»        0.00 
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TABLE    9 


SAMPLE 

OUTPUT 

-FROM    " 

•ROCKET" 

,    CASE 

II 

TIM 

<^.- 

FD3: 

'?^ii  VGl 

:-  nEIarT 

»33" 

I  ha  *• 

DJU  , 

A" 

!  1  '  A  I 

7.205 

5S5.14 

-555.13 

379.50 

33.81 

0.00 

35.12 

3.851 

coooc 

1  .i.i.J 

J  /  <L  •  o3 

-545.45 

7543.90 

^O'l    9(3 

33.04 

0.00 

34.33 

0.  0%j  t 

0. 0000 

7.246 

559.73 

-634.70 

7554.42 

380.95 

32.27 

0.00 

JO.  JJ 

8.851 

0.0000 

7.266 

546.35 

-524.11 

7564.68 

381.50 

0.00 

32.78 

8. 851 

O.OCCO 

7.287 

534.18 

-613.15 

7574.71 

382.24 

30.78 

0.00 

32.03 

8.851 

0.0000 

7. 307 

521.73 

-6C2.31    . 

7584.50 

352.86 

30.06 

i;.00 

8.851 

0.0000 

7.327 

509.51 

-591.39 

7594.07 

383.47 

29.35 

0.00 

jv>  56 

3.851 

0.0000 

■?    ni  0 
/.OtO 

497.51 

-580.60 

7603.40 

384.07 

28.65 

o.oo 

29.84 

S.851 

0.0000 

7.366 

-563.35 

wS7.DJ 

27.97 

-0.00 

29.14 

3.851 

0.0000 

7.389 

-559.22 

7621.41 

335.21 

27.30 

o.oo 

'2S.45 

,  3.851 

o.oooo 

7.409 

462.32 

-543.69 

7630.09 

335.76 

iO.&4 

0.00 

27.77 

2     rc> 

Q.OJl 

o.oooo 

7.430 

451.69 

-533.30 

7633.56 

0.00 

27.11 

6.351 

0.0000 

7.450 

440.77 

-525.02 

7546.83 

336.83 

25.36 

0.00 

ij.4c 

3.351 

430.06 

-517.90 

'£54.89 

337.44 

i7.  /'7 

o.cc 

25.32 

•j.SJi 

0.0000- 

7.^91 

419.55 

7652.76 

337.84 

0.00 

';C       '  Q 

3.35i 

'J  %  \}  -J  V  V 

■•CIO 

-49S.04 

7570.43 

,■     -            ;   •, 

0.  vU 

24.53 

3.351 

o.oooo 

1  .  J.jfc 

399.16 

-4cB.26. 

7677.31 

0(3w.  0^ 

o.-oo 

23.97 

3.351 

^     .^  A  A  A 

/•J/ J 

379.55 

-469.21 

^1.  C. 

coo 

22.30 

O.SJ. 

A      AAAA 

7,614 

3oG.71 

7708.01 

390.59 

20.72 

0.00 

21.53- 

-    nc( 

O.OJl 

0.0000 

7.655 

342.62 

-432.72 

7T4  a: 

391.42 

19.68 

0.00 

20.60 

W  •   U  w  « 

V.  1^000 

7.636 

'T-'C     1C 

392.21 

13.57 

0.00 

19.56 

3.851 

0.0000 

7.737 

305.59 

7,'Tj.Oo 

332.95 

11    71 

*  / .  /  * 

0.00 

i  J.  JO 

C.  O^i 

cooco 

7.778 

232.62 

-331.89 

7754.15 

393.55 

15.73 

0.00 

17.61 

n    nr  - 
3.  OJl 

A  AAA  A 
V.UVVV 

7.613 

277.30 

-365.95 

7764.68 

394.33 

15.90 

0.00 

16.53 

3.351 

0.0000 

7.860 

262.63 

-350.52 

7774.64 

394.96 

15.05 

0.00 

15.31 

3.351 

0.0000 

7.901 

248.58 

-335.56 

7704.07 

335.56 

14.24 

0.00  . 

.7.3; 

2     P^.t 

0.0000 

7.942 

ta  J  J.     A«J 

-321.03 

7792.99 

356.13 

13.47 

1 4 .  i  (■ 

^    0-  ' 

-V     ,-.  ^.  A  -. 

v.  W  v'Vv 

7.933 

^.>..    ^/ 

-307.06 

,'  :vl .  ^i 

-9i.b6 

'   -      TA 

0.00   ■ 

13.40 

w.  Qui 

.-.      A,:-  A  A 

3.024 

209.97 

-i!'j3.j0 

730:. 38 

397.17 

0.00 

12.56 

3.351 

0.0000 

3.055 

198.22 

-^oU. JO 

7316.90 

397.65 

11.34 

o.oo 

11.95 

5.351 

0.0000 

3.106 

186.99 

-267.70 

7823.99 

"^^  If) 

1^.69 

0.00 

i  1  ..iO 

3.351 

0.0000 

8.147 

176.23 

-255.44 

7330.57 

393.52 

10.07 

'  o.oo 

10.64 

3.S51 

0.0000 

8.188 

165.21. 

-221.51 

7E36.97 

398.92 

S.54^ 

0.00 

10.04 

8.351 

0.0000 

3.249 

155.54 

-133.82 

7S45.95 

-    399.50 

•■»  9.09 

coo 

9.40 

6.851 

o.oooo 

8.341 

145.01 

-98.55 

7351.72 

400.31 

*3.55 

0.00 

8.75 

8.851 

0.0000 

8.464 

133.96 

-84.15 

7874.50 

401.32- 

'7.91 

0.00 

£.10 

8.851 

0.0000 

S.62B 

120.86 

-76.13 

7893.98 

402.55 

7.14 

0.00 

7.31 

8.851 

0.0000 

8.792 

108.87 

-69.82 

7911.45 

403.66 

6.43 

0.00 

£.53 

8.851 

c.cooo 

8.956 

97.89 

-64.30 

7927.17 

404.65 

5.78 

0.00 

5.93 

6.651 

0.0000 

S.I19 

87.79 

-53.96 

7941.29 

405.56 

5.13 

coo 

5.32 

S.351 

0.0000 

5.324 

/O.  JD 

-52.70 

7956.89 

406.55 

4.51 

0.00 

4.63 

8.851 

0.0000 

9.529 

56.16 

-47.05 

7970.42 

407.41 

3.91 

coo 

3.851 

0.0000 

9.775 

55.39 

-40.86 

7934.27 

408.30 

3.27 

0.00 

3.37 

d.oji 

«'^       ■■  'i  r-  '. 

10.062 

44.65 

.  -34.30 

7qC7    cc 

409.14 

2.64 

0.00 

..72 

3.851 

.;.oooo 

10.389 

34.48 

-27.91 

n  A  /■»  0    c  *■> 

409.90 

.    2.04 

0.00 

2.10 

8.351 

0.0000 

10.799 

24.48 

0  V  i.  V  •  o  g 

410.61 

0.00 

:.50 

S.851 

I'.  GO  00 

11.311 

15.62 

-■:  •    -ii. 

5030.0? 

411.20 

0.91 

0.00 

v..  94 

3.S51 

•;       A    ■   A  A 
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TABLE    9     (CONTINUED) 
SAMPLE    OUTPUT    FROM    "ROCKET",    CASE    II 


AVERAGE  PRE33:J-£  FRCr 
HAJtlrui?!  rRESS'jSi  liAS 


1017.2  PSI 


O.OOi 


Q- 


w  a      Rflu 


u'.u 


■3i 


!  :;«£ 


PC 


7.573 

373.55 

J.JO 

7.5:4 

360.71 

3.30 

7.655 

1 ,'  •-.  '  1 

0.J4 

7.696 

7.737 

303.59 

3.41 

7.778 

292.62 

3.45 

T.ai9 

277.30 

j."0 

7.360 

262.63 

•^  c  ■ 

7.901 

243.58 

3.  jb 

7.942 

"    /"A 

7.932 

-■'-  nn 

3.64 

Z    ■•'•■  -^ .' 

209.97 

3.55 

S.065 

198.22 

0.  /i. 

a.iOb 

153.99 

2.77 

2.147 

176.26 

J.  31 

3. 1S3 
3.249 

166.21 

3.66 
5.91 

S.341 

1  iS    iM 

3.96 

3.464 

4.03 

2  u'/D 

12v.S& 

<  * 

t.  Li 

3.792 

;'>.  3.0/ 

4. 20 

3.956 

:7.5: 

4,29 

Q  MQ 

J  •  *  *  rf 

57.79 

-».oo 

9.324 

/a. Ob 

4.51 

9.529' 

6o.  lo 

4.54 

9.775 

4.31 

iU. vb2 

44.65 

5.02 

10.389 

34.48 

5.23 

10.799 

24.48 

5.56 

11.311 

15.52 

6.19 

4368.49 
4111.61 
3373.89 
3646.54 

w.»..  Of 


204.4.3c 

*  tf  • w ■ w^ 

<7SC  "^T 
.  /  C  C .  .J  / 

lu4d. 13 

I  -T  .-t  A  r,4 
iOV'J.rf4 

*  i w*..  Jv 

»  A  <  ^   Ql 

tV:; . 3 J 

395.70 

■>  -  ■  c 

650.20 

JJv. /O 


.17 
">"/-    Ifl 

224.00 

126.70 

43.51 


INT  PC 

6506.01 
6521.17 

bjju. JO 
£549.25 

Uuw...  ^^ 

•3^7i  5^ 


j'.'«.J.  ■•/ 

OJOC. ^^ 

2334. 4v 

G J  ^ / . ^u 

2654.13 

Obv/ . /j 

2432.33  ' 

S617.65 

2323.52 

*•  ^  *  T    A  .-I 

2153.17 

ofcjj.e/ 

CO>»t.iO 

6659.56 

<■    i*  f   r        C'^ 

6o£i.2; 

5690.31 

b/.;B.j.. 


6754.08 

6796.10 
6310.70 

DO..J.  Cw 

6539.96 
6552.94 
6365.02 
5375.28 


INT  F 

85251.34 
85425.01 
6oj83. j5 
56742.58 
36337.49 

C  /  V ,.  J  .  /  'J 

d  <  i  J  ,  •  C  'J 

J  /  ^ .' .  •  J  / 
o>  OwJ.^a 
o'toj. ;/ 
37537.70 

h/ wUV. «v 

'3/  'bb.  /O 

37347.75 

«^^/^«  ^^ 
U  <  w^u  a  U/ 

ao'.-jj.  w4 

«  «  A  i^  A   i  c 

90<.ww.  Tu 

rn'i,--  AC 
C2J:3.'JJ 

53593.0- 
c377£.54 
53932.61 
3907 1 . 26 

89464.61 
39574.05 
8965^15 
33731.01 
89779158 


:I5HT 


ur 


cSr 


IS? 


389.72 

1.300 

390.59 

1.28S 

jii.t.. 

4  ..jC 

3':  3. 66 

1  Oii 

''.r,.  1^- 

.  ^«.» 

w?-. Jj 

1  .^J^ 

39^.96 

1.219 

395.55 

i    ^  .*.  ^ 

■70:  •■: 

OJiJ.  ij 

396.66 

1.184 

397.17 

1.175 

?C7  C? 

1.166 

393.10 

1.156 

393:52 

1.146 

399:5$ 

■   'a 

M25 

.  4  .^ 

TtU.Ji 

1.097 

jl  i*n  cc 
"l....  Jw 

A-r7 

>»•  J.20 

...  ,  V  M  W 

.,  JO" 

405. 5S 

.   A  1  .', 

tOt. jj 

0.977 

.407.41 

0.94-1 

403.30 

0.392 

409,14 

0.325 

409,90 

C.734 

410.61 

'v.  JO  J 

411.20 

0.315 

332.90  221.76 
3S9.81  221.71 
39'J.67  221.64 

10 1  CA  '?•■>•  ^7 


''92 

;■■-*' 

221.42 

7-"5 

221.25 

w'  . 

221.2s 

www 

'?"- 

-■  ■-.  •  <  c 

396. 

13 

221.0s 

w'jw« 

"'• 

..  - .  -^  ^. 

397. 

O'*^ 

^-^ 

297. 
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i^^O.  B: 

7Z' 

14 

fc*.v  .  JT 

4(n.. 

5e 
lb 

220.58 

401. 

220.43 

402. 

Ov 

*.iv.iJ 

40  j. 

**— 

A'.-   AT 
Li.\>t  'J  J 

219.53 

i .-  - 

■T  v  ».  • 

'^ '  -  "" '" 

4:;c. 

4v 

219.53 

TV;  ,iJ  i.i3,  jO 

iC9,07  215.97 

409.36  218.77 

410.50  216, j5 

411.2.  218.33 
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TABLE  10 


SUMMARY  OF  RESULTS  FOR  THREE  GRAIN  DESIGNS 
D  =  16  IN,  L  =  40  IN 


GRAIN  3 


II 


III 


No. 

of    pts. 

web. 

in 

Required    r    (in/sec) 

P.F. 

Port 

Area,    in^ 

RRl, 

in 

<t>2 

<j>4 

SRI, 

in 

SR2, 

in 

^max 

,    psia 

Fmaxlbf 

F  Ibf 

Irji,  lb£-sec 

^burn  '  ^®^ 
total 


6 

8 

8 

2.5 

2.5 

2.5 

0.385 

0.385 

0.385 

2.218 

2.553 

2.314 

42.48 

42.12 

46.88 

2.750 

3.0 

3.25 

26.81 

26.81 

26.81 

10.03 

3.90 

4.54 

0.5 

0.25 

0.25 

0.25 
990. 

12,464. 
6,956. 

91,151. 
13.1 


0.25 

1,017. 

13,702. 

7,945. 

89,780. 

11.3 


0.25 
952. 

12,749. 
7,738. 

87,690. 
11.2 


65 


3      4   5 

TIME  ,  sec: 


6 
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FIGURE  7. 


PRESSURE-TIME  PROFILES  FROM  "ROCKET' 
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X.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

An  interactive  computer  code  for  the  preliminary  design 
of  solid  rocket  propellant  rocket  motors  ("SPRMD")  was  suc- 
cessfully developed  and  its  use  was  demonstrated  through  a 
design  example.  The  code  incorporated  several  existing  PC 
codes  and  AGARD  recommended  empirical  expressions  for  losses 
in  aluminized  propellant  systems. 

B.  RECOMMENDATIONS 

Further  modifications  to  improve  the  capabilities  of  the 
code  are  essential  in  order  to  make  "SPRMD"  a  more  effective 
preliminary  design  tool. 

-  There  is  only  one  grain  configuration  (radially  spoked 
star  grains)  incorporated  into  "SPRMD." 

It  is  necessary  to  add  other  options  for  grain  design 
in  order  to  increase  the  grain  design  flexibility.  One 
possibility  would  be  to  incorporate  the  Grain  Design 
and  Internal  Ballistics  (GDB)  module  of  Solid  Per- 
formance Prediction  Code  (SPP)  (Ref.  3). 

-  "SPRMD"  also  does  not  currently  incorporate  any  plot- 
ting capabilities.  There  is  definitely  a  need  for  this 
program  to  have  plotting  capabilities  for  the  grain 
design  results  and  the  thrust  and  pressure  versus 
time  profiles. 
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"FLYIT"  was  not  incorporated  into  the  program,  although 
"SPRMD"  was  structured  to  have  flight  simulation  capa- 
bilities. It  is  also  recommended  that  "FLYIT"  be 
incorporated  into  "SPRPID"  as  another  enhancement. 
"SPRr4D"  is  a  primative  type  of  preliminary  solid  pro- 
pellant  rocket  motor  design  code.  Further  enhance- 
ments, such  as  propellant  searching,  screen  editing, 
and  increased  solution  speed  for  "MICROPEP"  and 
"ROCKET"  would  make  this  code  more  user  friendly. 
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program  SPRMD- 

c  * 

c  SOLID  PROPELLANT  ROCKET  MOTOR  DESIGN  * 

C  THIS  PROGRAM  (SPRMD)  SOLVE  THE  PRELIMINARY  DESIGN  PROBLEM* 
C  OF  SOLID  PROPELLANT  ROCKET  MOTOR.  SUBROUTINES  INCLUDES:  * 
C  (1)  INPUTO:  BASIC  INPUT  FOR  DESIGN  * 

c  (2)  INPUTl:  CHEMICAL  COMPOSITION  INPUT  * 

C  (3)  MICROPEP:  PERFORME  THE  PROPELLANT  EVALUTION  * 

C  (4)  PFE:  SIZE  THE  THROAT  OF  NOZZLE  PORT  AREA  OF  GRAIN  ETC* 
c  (5)  LOSCF:  CALCULATE  THE  THRUST  COEFFICIENT  * 

C  (6)  INPUT2:  STAR  GRAIN  INPUT  * 

C  (7)  GRAINS:  PERFORM  THE  SPOKED  STAR  GRAINS  DESIGN         * 
c  (8)  RCKTIN:  INPUT  THE  DATA  FOR  THRUST  AND  PRESSURE  PRO-   *  • 
c  FILE  CALCULATION  * 

c  (9)  ROCKET:  COMPUTED  THE  THRUST  AND  PRESSURE  PROFILE  * 
c*********************************************************** 
c 

real  f avg , tb , 1 , d , meop , pc , rrqd , pamb , tamb , h , Id , webrqd , wf , m 
common  /micrp/  a(12, 12) ,kr(20) ,amat (10, 12) , jat (12) , in,is, 

lfie(10,6) ,ie(10,6) ,alp(12) ,w27,n,dh(10) ,rho(10) , 

2wate(10) ,wl(6) ,w43 , ig,np,vnt (201) ,w47,name,ser, 

3flobr,itag(100) ,wing(10) 
c     integer 

character*10  yes 
c« ********************************************************** 

c  * 

c  BASIC  INPUT  FOR  MOTOR  DESIGN  * 

c  * 

c* ************* ********************************************* 

1000  call  inputO ( f avg , tb , 1 , d , meop , pc , rrqd , pamb , taunb , h , f min , tpmax , tpnom, 
*tpmin,pik) 
write(*,l) favg 

1  format (ix, 'AVERAGE  THRUST,  LBF ',£12.1) 

write ( * , 2 ) tb 

2  format(lx, 'BURN' TIME,  SEC ',£12. 2) 

write (*, 12) fmin 

12  format (ix, 'MINIMUM  THRUST  REQUIRED,  LBF ',£12. 6) 

write(*,3)l 

3  format (ix, 'MOTOR  LENGTH,  INCH ',£12. 2) 

write(*,4)d 

4  format  (ix, 'MOTOR  DIAMETER,  INCH ',fl2.2) 

write (*, 5) meop 

5  format(lx, 'MAXIMUM  EXPECTED  OPERATING  PRESSURE...  ',fl2.1) 
write  (*,  13)  tpmajt 

13  format  (Ix, 'TPMAX,  F i . . . .' ',£12.1) 

write(*, 131) tpnom  * 

131  format  (Ix, 'TPNOM,  F.... •, *  ,tl2.1) 

write(*,132)tpmin 

132  format ( ix, 'TPMIN,  F ',fl2.1) 

write (*, 14 )pik 

14  format(lx, 'TEMPERATURE  SENSITIVITY,  %/F ',£12. 3) 

write(*, 6)pc 

6  format (ix, 'NOMINAL  CHAMBER  PRESSURE,  PSI •,fl2.1) 

write (*, 7) rrqd 

7  format(lx, 'REQUIRED  BURNING  RATE,  IN/SEC ',fl2.3) 

write(*,8)h 

8  format (ix, 'ALTITUDE,  FT ',fl2.0) 

write (*, 9) pamb 

9  format (ix, 'AMBIENT  PRESSURE,  PSI ',£12. 2) 
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write(*,ll) tamb 

.   11   format (ix, 'AMBIENT  TEMPERATURE,  R ',fl2.1) 

write(*,10) 
10    format (Ix, 'ANY  CORRECTION  OF  BASIC  INPUT?  Y-YES,  N-NO...  ') 

read (*, 20) yes 
.20    format (alO) 

if (yes.eq. 'y' .or.yes.eq. 'Y')  go  to  1000 
c 
. c**************** ************************************* ****** 

c  '  *.         • 

c  INPUT  THE  DATA  NEEDED  BY  MICROPEP  * 

c  * 

c*********************************************************** 

c 

c  •  •      . 

1010  call  inputl(pc,pamb) 

write (*, 40) 
40    format (ix, 'Are  there  any  corrections  of  input  for  micropep?,y-yes 
*n-no  • ) 

read (*, 20) yes 

if (yes.eq. 'Y' .or.yes.eq. 'y* )  go  to  1010 
c 

c** ***************************************  ****************** 
c  * 

c  COMPUTE  THERMAL  EQUILIBRIUM,  PROPELLANT  CHARACTERISTIC,  * 
c  AND  ETC. (USING  MICROPEP,  UNTILL  GOOD  RESULTS  ARE  REACHED)* 
C  * 

c****** *********************************************** ****** 

c 

call  micropep  . 
write(*,50) 
50    format(lx, 'Are  the  results  good?  y-yes  n-no   •) 
read (*, 20) yes 
if (yes.eq. 'y' .or.yes.eq. 'Y')  then 

go  to  1020 
else 
55       write (*, 60) 
60       format  (Ix,  '  (1)   Change  the  Basic  Inputs  of  motor.')- 

write(*,70) 
70       format (Ix, ' (2)   Change  the  Input  data  for  Micropep.') 

write(*,80) 
80       format(lx, '     type  1,  or  2,  if  1  is  chosen  it  is  necessary') 

write(*,90) 
90       format(lx,'     that  Micropep  shall  be  run  once  again! 1 1!') 

read (*, 100) i 
100      format (il) 

if (i.ne. l.and.i.ne.2)   then  '  ' 

write(*,*)'      Error  inputi^l !  1,  or  2  must  be  typed' 
go  to  55  - 
else 

if(i.eq.l)  go  to  1000 
if (l.eq.2)  go  to  1010 
endif 
endif 
c* ********************************************************** 
c  * 

c  CALCULATE  THE  LOSS  PERFORMANCE  &  COMPUTE  THE  DATA  FOR  * 
C  GRAIN  DESIGN  * 

C  * 

c*********************************************************** 
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c 

1020  call  pfe(favg,tb,l,d,pc,cfc,ath,etacf ,m,vp,vl,dpe,aport, 
*pf ac , isp , cf th , webrqd , wf , tc , rrqd , ae , etac , xx) 
write (*, 50) 
read(*,20)yes 
if (yes.eq. 'y' .or.yes.eq. 'Y')  then 

go  to  1030 
else    • 
75    write (*, 60) 
write(*,70) 
write(*,71) 

71  fonnat(lx, ' (3)   Change  nozzle  data.') 
write(*,72)  , 

72  format ( Ix, '     type  1,  2,  or  3   •) 
read(*,100)i 

if (i.lt.l.and.i.gt.3)   then 

write (*,'>)•  Error  input!  1  1,  2  or  3  must  be  typed' 

go  to  75 
else 

if(i.eq.l)  go  to  1000 

if(i.eq.2)  go  to  1010 

if(i.eq.3)  go  to  1020 
endif 
endif 
c* ********************************************************** 

c  * 

c  CHOOSE  THE  APPROPRIATE  GRAIN  DESIGN  * 

c  NOTE:  OPTION  2,  3,  4,  5,  6,  7  ARE  NOT  AVAILABLE.  INPUT  * 

c  YOUR  OWN  DESIGN  BY  SUBSTITUTE  THE  WRITE  AND  GOTO  STATE-  * 

C  MENTS  IN  THE  CORRESPONDING  BLOCK  * 

C  * 
c* ********************************************************** 

1030  write(*,110) 

110   format (Ix, 'Choose  one  of  the  following  grain  configurations') 

write(*,120) 
120   format ( Ix, ' (1)  Star  grain') 

write(*,130) 
130   format ( Ix, ' (2)  Optional  grain  •) 

write(*,140) 
140   format(lx, ' (3)  Optional  grain  •) 

write (*, 150) 
150   format ( Ix, ' (4)  Optional  grain  •) 

write (*, 160) 
160   format ( ix,  '  (5)  Optional  grain  •)  '' 

write (*, 170)  .  , 

170   format(lx, ' (6)  Optional  grain  ]) 

write (*, 180) 
180   format (Ix, • (7)  Optional  grain  ') 
write(*,190) 

190   f ormat ( Ix , ' Type  1,  2,  3,  etc.') 

read(*,*)i 

if (i.lt.l.or.i.gt.7)  then 
write(*,200) 

200     format (Ix, 'Error  input 1 !  must  type  1,  2,  3, ,7') 

go  to  1030 
else 
if(i.eq.l)  then 
rr5=d/2.0 

write (*,*) 'Run  "DESIGN  MODE"   • 
1050    call  input2 (pfac,aport, webrqd, rr5,l,tb, XX) 
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call  grains 

write(*,*) 'Are  the  results  are  acceptable?  y-acceptable  n-not 
♦acceptable  type  y  or  n. . ' 
read(*, 20) yes 
if (yes.eq. 'n' .or.yes.eq. 'N')    then 

go  to  1050 
else 

go  to  1055 
endif 
1055    rr5-d/2.0 

write(*,*) 'Run  the  "BURNBACK  MODE"  to  get  web  vs.  burning  area' 
call  input2 (pfac,aport,webrqd,rr5,l, tb,xx) 
call  grains 

write(*, *) 'Are  the  results  are  acceptable?  y-acceptable  n-not 
♦acceptable  type  y  or  n. . ' 
read(*,20)yes 
if (yes.eq. 'n' .or.yes.eq. 'N' )  then 

write(*,*) 'Do  you  want  to  >:ry  another  star  grain?  y/n   ' 

read (*, 20) yes 

if (yes.eq. 'n' .or.yes.eq. 'N')  then 

go  to  1030 
else 

go  to  1050 
endif 
else 

go  to  1060 
endif 
else 

if(i.eq.2)  then 
write(*,280) 
280         f ormat ( Ix ,' Input  your  own  grain  design  package') 
go  to  1030 
else 

if(i.eq.3)  then 
write(*,290) 
290  format (Ix, 'Input  your  own  grain  design  package') 

go  to  1030 
else 

if(i.eq.4)  then 
write(*,300) 
300  format (Ix, 'Input  your  own  grain  design  package') 

goto  1030 
else  '    .... 

if(i.eq.5)  then 
write(*,310) 
310  f ormat ( Ix ,' Input  your'own  grain  design  package') 

go  to  1030  »" 

else  ' 

if(i.eq.6)  then 
write(*,320) 
320  format (Ix, 'Input  your  own  grain  design  package') 

go  to  103  0 
else 

if(l.eq.7)  then 
write(*,330) 
330  format (Ix, 'Input  your  own  grain  design  package') 

go  to  103  0  i  .. 

endif 
endif 
endif 
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endlf 
endif 
endif 
endif 
endif 
c 

c  .  * 

C  INPUT  THE  DATA  FOR  THRUST  &  PRESSURE   VERSUS  TIME  PROFILE* 
c  CALCULATION  * 

c  * 

c*********************************************************** 
1060  write(*,*) 'The  following  processes  is  to  calculate  the  "THRUST  VS, 
*  TIME"  and  "PRESSURE  VS.  TIME"  profiles'  • 

va- (d/2 . 0) **2 . 0*3 . 14159*1 
rrqd»webrqd/tb 
write (*, 106) rrqd 
106  format (Ix, 'Required  burning  rate  for  new  web,  IN/SEC.  ',f6.4) 
1070  call  rclctin(ath,ae,va,tamb,painb,pc,etacf ,rrqd,pik,etac) 
write (*,*) 'Any  change  for  the  input?  y-yes,  n-no    ' 
read (*, 20) yes 
if (yes.eq. 'y' .or.yes.eq. 'Y')  then 

go  to  1070 
else 

go  to  1080 
endif 
c 

c**** ******************************************** *********** 
c  * 

c  COMPUTE  THE  PRESSURE  VS.  TIME   &   THRUST  VS. TIME  PROFILE  * 
!       c  OF  ROCKET  MOTOR  * 

"       C  .* 

c************************ ************************** ********* 
c 

1080  call  rocket 
c 
!    '       write(*,340) 

340   format (Ix, 'Are  the  results  acceptable?  y-yse,  n-no •) 

read (*, 20) yes 

if (yes.eq. 'y' .or.yes.eq. 'Y') then 
I  go  to  9999 

I  endif 

i  write(*,350) 

i      350   format (Ix, 'Do  you  want  to  proceed?  y-yes,  n-no ') 

!  read (*, 20) yes 

I  if (yes.eq. 'N' .or.yes.eq. 'n')  tKen 

write  (*,*)•  Terminated  on  the*  request  of  user  ' 
go  to  9999 
endif 
1040  write(*,360) 
360   format (Ix, 'Suggest  to  change  the  following  data') 

write(*,380) 
380   format (Ix, ' (1)  Change  the  Grain  Design.') 

write(*,390) 
390   format (Ix, ' (2)  Change  the  Ingredients  of  the  propellant. •) 
■  .  write(*,400) 

400   format (ix, ' (3)  Change  the  Basic  Design  Requirements.') 
write (*,*)' (4)  Re-run  the  "ROCKET"' 
write (*,*)'  Type  1,  2  or  3  for  selection' 
read(*,100)i 
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If (l.ne.l.and.l.ne.2.and. i.ne.3)  then 

write (*,*)'  Error  input  you  have  to  type  1,  2  or  3 ' 
go  to  1040 
else 

if(i.eq.l)  then 

go  to  1030 
else 

if(i.eq.2)  then 

go  to  1010 
else 
if(i.eq.3)  then 

go  to  1000  I 

else 

if(i.eq.4)  go  to  1070 
endif 
endif 
^ndif 
endif 
c**** **************************  ************  ***************** 
c  * 

c  INPUT  THE  TRAJECTORY  DATA   (FUTURE  WORK)  * 

c  * 

c* ****** **************************************************** 

c 

c     call  flyin 

c 

C* *************************************************'*  ******** 

c  * 

C  CALCULATE  THE  TRAJECTORY  OF  ROCKET  MOTOR  (FUTURE  WORK)    * 
c  * 

c*********************************************************** 

c 

c     call  flyit 

9999  end 


:« 

* 
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subroutine   Input 0 (favg, tb,l,d,meop,pc,r,paaib, tamb,h, f min , tpmax , 
* tpnom , tpmin , pik) 

c  * 

c  THIS  SUBROUTINE  IS  FUNCTION  AS  A  BASIC  INPUT  FOR  SOLID    * 
C  PROPELLANT  ROCKET  MOTOR  DESIGN  * 

C  .  * 

c* ******************************  **************************** 

c 

c  varizdale  definition; 

e     *********  (ballistic  performance)  ********* 
c     favg    :  average  thrust  (F12.6) 
c     tb      :  burning  tine  of  motor  (F12.6) 
c     fmin    :  minimvim  thrust  (F12..6)  * 

c     *********  (mission/vehicle)  *************** 
c     1       :  motor  length  (F12.6) 
c     d       :  motor  diameter  (F12.6) 
c     meop    :  maximum  expeted  operating  p^^'ssure  (F12.6) 
c     va      :  motor  volume  avaiable  for  propellant  (F12.6) 
c     tpmax,  tpnom,  tpmin 
i      c     pmaxtp  :  meop  *  (1-reproducibility  tolerance)  (F12.6) 
i      c     pc      :  norminal  chamber  pressure  (F12.6) 
I      c     rptol   :  reproducibility  tolerance  (F12.6) 
;      c   ■  r       :  burning  rate  (F12.6) 
c     web     :  web  of  grain  (F12.6) 

c     pik     :  temperature  sensitivity  of  propellant (F12. 6) 
c     h       :  altitude (F12. 6) 

real  favg , tb , fmin , 1 , d , meop , va , dltp , pmaxtp , pc , rptol , r , web , h 
c     open  a  file  which  can  hold  the  basic  input  data 

open ( 2 , f ile= ' inputO . dat ' , access= ' sequential • , status- ' unknown • ) 
c     input  basic  data 

write(*,*) 'All  input  data  in  decimal  format' 

write (*,*) 'Specify  the  average  thrust,  LBF • 

read  (*,10)favg 
10    format (f 12. 6) 

write(2,10) favg 

write (*,*) 'Input  burning  time  of  motor,  SEC • 

read(*,10)tb 

write (2, 10) tb 

write (*,*)' Input  minimum  thrust  required,  LB ' 

read (*, 10) fmin 

write (*,*) 'Input  motor  length,  INCHES • 

read(*,10)l 

write (2, 10)1 

write (*,*) 'Input  motor  diameter,  INCHES • 

read(*,10)d  a 

write (2, 10) d  * 

write (*,*) 'Input  maximum  expected > • 

write (*,*) 'operating  pressure,  PSI ' 

read  (*,  10)  meop  '* 

write (2, 10) meop 

I  write (*,*) 'Input  Tpmax,  F • 

♦   read (*, 10) tpmax 
i  write(2, 10) tpmax 

write (*,*) 'Input  Tpnom,  F • 

read ( * , 10 ) tpnom 

write ( 2 , 10 ) tpnom 

write(*,*) 'Input  Tpmin,  F. ' 

read(*, 10) tpmin 

write ( 2 , 10 ) tpmin 


75 


APPENDIX 


write(2,10)dltp 

write(*,*) 'Input  reproducibility  tolerance,    % ' 

read ( * , 10 ) rptol 

write ( 2 , 10) rptol 

rptol=rptol/ 100.0 

pinaxtp=meop*  ( 1 . 0-rptol ) 

rptol=rptol*100.0 

r=(d/4)/tb 

write (*, 20) r 

write(2,10)r 

20.    format (Ix, 'Estimated  repuired  bum  rate ',fl2.3, 

*'  INCH/SEC) 

write(*,30) 
c     The  following  calculation  of  ambient  temperature  and  pr&ssure  at 
c     certain  altitude  is  directly  copied  from  NWC  Technical  Memorandum 
c     4757  subroutine  atmstd  based  on  U.S  Standard  Atmosphere, . 1976. 
30    format (Ix, 'Input  the  required  design  altitude  of  rocket  motor  in  F 
*T  

*•) 
read(*,10)h 
write(2,10)h 
if (h-36152.)l,l,2 

1  tamb=518.67-0.003559969*h 

ro=0.00237696*exp(-(0. 028800339* (h*. 001) +0. 000132317* (h*. 001) **2)) 
go  to  7 

2  if (h-65825.)3,3,4 

3  tamb=389.97 

ro=0. 00237696*0. 29697*exp(-0. 04771923* (h-36152.)*. 001) 
go  to  7 

4  if (h-105511.)5,5,6 

5  tamb=354. 142+0. 0005442926*h 

ro=0. 00237696*0. 0718594*exp(-(0. 049130937* (h-64825.)*. 001 
&-0. 000034655* ((h-6S825.)*. 001) **2)) 
go  to  7 

6  tamb=251. 4114+0. 00151794*h 

ro=0. 00237696*0. 0107993 *exp(-(h-105511.)/21320.) 

7  pamb=11.9203*ro*tamb 
write(2,10)pamb 
write(2,10) tamb 

c     end  of  calculation  of  tempressure  and  temperature  at  design  alt. 
c 

write(*,*) 'Based  on  available  requirements,  make  an  initial' 

write (*,*) 'selection  of  propellant.  Then  input  the  temperature' 

write (*^*) 'sensitivity  of  the  propellant,  %/F '  • 

read<*,40)pik 

write{2,40)pik  ? 

40    format(fl2.6)  * 

pik=pik/100.0 

pc=pmaxtp*exp (pik* (tpnom-tpmax) ) 

pik=pik*100. 0 

write(2, 10) pc 

close (2 , status= ' keep ' ) 

return 

end 
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i  subroutine   inputl(pc,parab) 

;         =•  * 

■         C  THE  PRIMARY  PURPOSE  OF  THIS  SUBROUTINE  IS  TO  ESTABLISH  * 

c  The  DATA  FILE  WHICH  CAN  BE  ACCESSED  BY  MICROPEP  SUB-  * 

C  TINE  * 

I         c  •  *        ■   . 

c* ******************************************** ************** 

i  ■  c 

1  •  common  /micrp/  a(12, 12) ,kr(20) ,anat(10,12) , jat(12) ,in,is, fie(10,6) 

I  l,ie(10,6)  ,alp(l2)  ,w27,n,clh(10)  ,rhoaO)  ,wate(10)  ,wl(6)  ,w43,ig,np, 

2vnt(201) ,w47, name, ser, floor, itag(lOO) ,wing(10) 
\  common/chara/block,je,aspec, specie 

!  character*8  specie (200)  «  • 

character*30  block(lO) ,of ile 
j  character*2  je(10,6) ,aspec(12) 

common/moon/tstest , te , irun , iounit , iend 

character*30  blok(lO) 

character*10  user 

dimension  jie(10,6) 

integer  nic,niu, irun, alt, kr 

real  denexp 
c  "      - 

c     open  the  file  as  formatted  sequential  as  unit  3 

open  ( 3 ,  f  ile=  •  input .  dat '  ,  access=» '  sequential  • ,  status^ '  unknown  • ) 
c  input  the  case  name  from  the  keyboard 

write(*,*)'  Start  to  input  the  selected  propellant  data     .' 

write(*,*)'  • 

write(*,l) 
1  format (Ix, 'Specify  the  case  name  (format  AlO) •) 

read(*, • (a) • )user 

write (3, 10) user 
10    format (lalO) 

irun=2 
c  input  the  number  of  ingredients,  supplied  ingredients,  and  runs 

write (*,*) 'Input  the  number  of  ingredients?  (format  15)....' 

read(*,*)nio 

write (*,*) 'Input  the  number  of  user  suppplied  ingredients  (15)...' 

read(*,*)niu 

wr it e ( 3 , 2 5 ) nic , niu , irun 
20    format(2i5) 
25    format(3i5) 

c  Input  density  exponent  which  is  used  in  conjunction  with  control 
c  option  3. 

write  ■(*#*) 'Input  the  density  exponent:' 

write (*,*) 'For  tactical  missiles  or  first  stages  use  1.0' 

write(*,*) 'Use  0.7  for  second  stagps  and  0.2  for  third  stages' 

write (*,*)' (format  FIO.O) ' 

read (*, 4 0 ) denexp 

wr i te ( 3 , 4 0 ) denexp 
40    format(fl0.5) 
c 

c  input  control  option  store  it  in  kr(20) 
C     print  out  the  options  on  the  screen 

write (*,' (a) ') 'There  are  five  options  listed  below  for  your' 

write(*,*) 'selection.  Type  1,  2,  3,  4,  5—' 

write (*,*) 'Normal  input  for  design  is  3' 

write(*, ' (a) ') '    1.  For  chamber  and  exit  calculations. ' 

write(*, ' (a) ' ) '   2.  For  chamber  only  calulations. ' 

write(*, ' (a) ' ) '    3.  For  boost  performance  and  nozzle  design.' 
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45 


write(*, • (a) 
write(*, • (a) 
write(*, • (a) 
write(*, • (a) 
write (*, • (a) 
write (*,*) 'Use 
read(*,45)alt 
format (il) 
if  (alt. It. 5) 
kr(2) 
kr(4) 
kr(5) 
kr(6) 
kr(8) 
if 


4.  For  equilibrivun  calculations  at  specified 
and  P.  Use  T  in  Kelvin  at  Pe  location  and 
in  psia  at  the  normal  Pc  location. • 

5.  Other  options  (  Check  README  for  special' 
alternatives) . • 

format  II  • 


then 
«  0 
-  0 
»  0 
=  0 
=  0 


T' 
P' 


(alt.eq.l) 
kr(l) 
kr(3) 
kr '7) 

else 

if  (alt.eq 


then 
-  0 
=  0 
=  0 


2)  then 

kr(l)  «  1 

kr(3)  -  0 

kr(7)  -  0 

else 

if  (alt 


else 


.eq.3) 
kr(l) 
kr(3) 
kr(7) 

kr(l) 
kr(3) 
kr(7) 


then 
-  0 
=  1 
=  0 

»  0 
"  0 

=  1 


50 


60 

c 
c 

70 


80 

90 

100 

120 

c 


endif 
endif 

endif 

write(3,50) (kr(i),i  =  1,8) 

format ( 8 il) 
else 

write(*, • (a) •) •  Inout  your  requirements  for  calculation — ' 
read(*,60) (kr(i) ,i'=  1,15) 
write(3,60) (kr(i),i  =  1,15) 
format  (15il) 


endif 

input  the  serial 


number  of  ingredients  found  in  JANNAF.DAT  table 


', i2 ,'• ingredients  needed  for  calculation') 


write(*,70)  nic 
format ( Ix, ' Input 
do  100  i  «  l,nic 

write(*,80)i 

forTnat(lx, 'Serial  #  of  number   •,i5,'    ingredient (15) .. .        •) 

read(*,90)itag(i) 

fonnat(i5) 
continue 

write  (3,120) (itag(i) ,  i  »  l,nic) 
format(10i5) 
Input  user  supplied  ingredients 

(niu.gt.O)  then 

170  i  =  l,niu 

wrice(*, ' (a) •) 'Specify  ingredient  name  use  caps  lock  (A30) 

read(*,130)blok(i) 


if 
do 
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write(*^' (a) •) 'Specify  composition  of  ingredient  max  =  6' 
write(*, ' (a) ') '  use  (6(13, A2))' 
write(*, '(a) •) •  ###**###**###**###**###**###**• 
read(*,140) (jie(i,l) , je(i,l) ,1-1,6) 

write(*, • (a) •) •  Specify  heat  of  formation,  and  density  (f5.0,f6 
*.0)     • 

read(*,15Q)dh(i) ,rho(i) 

write(3,160)blok(i) , (jie(i,l) ,je(i,l) ,1-1,6) ,dh(i) ,rho(i) 
130      format(a30) 
140      format(6(i3,a2) ) 
150      format(f5.0,f6.0) 

160      format(a30,6(i3,a2) ,f5.0,f6.0)  ^        . 

170   continue 
else 

write  (*,'(a)')'  No  user  supplied  ingredients.  • 
endif 
if  (irun.lt. 2)  then 

irun  =2 
endif 
status=0 
do  220  i  =  l,irun 

if (status. eq.O)  then 
wl(5)=250 
c       write  (*,'(a)')'  Specify  the  desired  chamber  pressure  (F12.4)  • 
c       read  (*,180)wl(5) 

wl(6)=pamb 
c      .  write  (*,'(a)')'  Specify  the  exhaust  pressure  (F12.4) • 
c        read  (*,190)wl(6) 
ni=nic+niu 
do  175  j  «l,ni 
write(*,172)itag(j) 
172      format (ix, 'Input  the  weight  %  of  ingredient' ,i5, • (FIO. 4) ...• ) 

read (*, 19 0 ) wing ( j ) 
175      continue 
180      fonnat(fl0.4) 
190      format(fl0.4) 

write(3,210)wl(5) ,wl(6) , (wing(j) ,j-l,10) 
210      format(12fl0.4) 

status-1  . 

else 

wl(5)=pc 

write(3,210)wl(5) ,wl(6) , (wing(j) , j-i;i0) 
endif  •  -'    . 

220   continue 

close  (3)  .^ 

return  •  \ 

end 
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subroutine  pf e ( f avg , tb , 1 , d , pc , cf c , ath , etacf , m , vp , vl , dpe , aport , 
+pf ac , isp , cf th , webrqd , wf , tc , r , ae , etac , xx) 

c  * 

C  THE  MAIN  PURPOSE  OF  THIS  SUBROUTINE  IS  TO  CALCULATE  THE   * 
c  PERFORMANCE  LOSS  AND  THE  DATA  FOR  GRAIN  DESIGN  * 

c  * 

c***********************.*  *********************************** 

c  variable  definitions: 


c  favg: 

c  cfth: 

c  cfc 

c  pc 

c  d 

c  1 

c  athth: 

c  ath 

c  atht 

c  dthth: 

c  etacf: 

c  dif 

c  va 

c  tb 

c  vp 

o  vl 

c  aport: 

c  dpe 

c  pwet 

c  j 

c  rhop 

c  mw 

c  mf 

c  webrqd: 

c  wf   : 

c  pfac  : 

c  character*80  fneune 

c 


Average  thrust  from  basic  input 

thrust  coefficient  from  micropep 

corrected  cf 

chamber  pressure  from  basic  input 

motor  diameter  from  basic  input 

motor  length  from  basic  input 

theoretical  throat  area 

throat  area 

tempory  result  of  throat 

theoretical  throat  diameter 

thrust  coefficient  efficiency 

for  convergence  test  of  throat  size 

vol  avaiable  for  propellant 

burn  time 

propellant  vol. 

volumetric  loading  factor 

port  area 

equivalent  circular  port  diameter 

wetted  perimeter 

throat-to-port:  area  ratio 

propellant  mass  density 

molecular  weight  of  propellant 

mole  fraction  of  condensed  phase 

required  web 

web  fraction 

perimeter  factor 


real  favg , cfth , cf c , pc , d , 1 , athth , ath , dthth , etacf , dif , tb , vp , vl , va , j , 
+aport, dpe, rhop, a, nn,rstd,rpc, web, webrqd, wf,m, pwet, pfac, isp, ispm,at 
+ht , mw , tc , ispth , Id , mf 
character*!  yes 
The  purpose  of  this  common  instiruction  is  to  transfer  w4 3  (propellant 
density)  from  MICROPEP:  SUBROUTINE  PUTIN (LE) 

comraQn  /micrp/  a(12,12) ,kr (20) ,amat(10,12) , jat(12) ,in,is, 
lfie(10,6) ,ie(10,6) ,alp(12) ,w27,n,dh(10) ,rho(10) , 
2wate(10) ,wl(6) ,w43, ig,np,vnt (201) ,w47, name, ser 
3,floor,itag(100) ,wing(10)       * 
The  purpose  of  this  common  instruction  is  to  transfer  the  CHAMBER  PRE- 
SSURE &  MOLECULAR  WEIGHT  from  MICROPEP  SUBROUTINE  OUT. 

common  /tcf/  tfc,rmwtc 
Transfer  SPECIFIC  IMPULSE,  CHARACTERISTISC  VEL. ,  CODFF.  OF  THRUST  & 
OPTIMUN  EXPANSION  RATIO  from  MICROPEP  SUBROUTIN  DESING. 

common  /socf/  spi (2) ,cf (2) ,oex(2) 
calculate  the  theoretical  throat  area 

pi»3. 14159 
transfer  cfth  from  "MICROPEP" 

cfth=cf (2) 
initialize  cfc 
cfc=cfth 
athth=favg/ (cfth*pc) 


fi 
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c  initialize  ath=athth 

ath=athth 
o  calculate  theoretical   throat  diameter 

dthth=sqrt (ath/pi) *2 

write (*, 888 )cfth,athth,dthth 
888   format(lx, 'Cfth=',f5.3, '   Athth-' ,f6.3, •   Dthth- • , f 5 . 3 ) 
c  calculate  the  etacf 

wr it e (*,*) 'Input  the  nozzle  half  angle' 

write (*,*) 'alpha,  DEGREES  (F5.2) ' 

read (*,*) alpha 

pi=3. 14159 

alpha= ( alpha/18 0. 0 ) *pi 

write(*,*) 'Input  the  nozzle  exit  half  angle' 

write(*,*)  'thetaex,  DEGREES  (F5.2) '   , 

read ( * , * ) thetaex 

thetaex=(thetaex/180.0) *pi 

write(*,*) 'Input  radial  erosion  rate  of  the  throat  in  IN/SEC 

write (*,*) 'assumed  to  start  at  time  zero  (FORMAT  O.XXXX)....' 

read(*, *) erosr 

write (*,*) 'Input  the  mole  fraction  of  condensed* 

write(*,*) 'phase  in  (MOLE/100  g) ' 

read(*,*)mf 

write (*,*) 'Input  the  material  of  nozzle:' 
•  write (*,*)• 1-steel  nozzle,  0-ordinary  nozzle ' 

read  (*,*)nozle  .  . 

atht=ath 

do  9  i=l,50 
C  calculate  the  loss  mechanism 

call  loscf ( etacf ; tb , ath , pc , alpha , thetaex , erosr , mf , nozle) 
c  correct  the  trust  coefficient 

cfc=cfth*etacf 
c  calculate  the  corrected  throat  area 

ath=f avg/ (cf c*pc) 

write(*, 8888) etacf ,cfc, ath 
8888   formatdx, 'Etacf  =',f  12.3,  '   Cfc=' ,  fl2.3,  '   Ath-',fl2.3)  .. 
c  calculate  the  convergence  athth  and  ath 

dif =abs ( (atht-ath) /atht) *100 . 0 

if  (dif. gt. 1.0)  then 
atht=ath 

else 

go  to  91 

endif 
9     continue  ^  ■       • 

91    va=pj,*(d/2)**2*l 

etac=0.93 

vol=pi*(d/4)**2*l 
c  Density  and  molecular  weight  com  froiij  "MICROPEP" 

rhop=w43        .  » 

mw=rmwtc 

tc=tfc 

ispth=spi(2) 

write (*, 98) rhop,mw,tc, ispth 

98  formate  Rhop=-' ,f  12.3,  '   MW-',fl2.3,'   Tc-',fl2.1,'   Ispth= ' ,  f  12 . 2 

*) 

99  isp=ispth*etacf*etac 
m=favg/isp 

vp=favg*tb/ (isp*rhop) 
vl=vp/va 

write (*, 87 )VP,VA,VL 
87  formate    VP=',F12.1,'      VA=',F12.1,'      VL-',F12.3) 
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dpe"sqrt ( 1 . O-vl ) *d 

aport=(pi/4)*d**2*(l-vl) 

tres=(aport*l)*pc*mw/(m*12*1545*(tc+459.6)) 
i  tres«tres*1000. 0 

write(*,88) tres 

I         88     formate  RESIDENCE  TIME  IN  MSEC •,fl2.3) 

;  write(*,89)etac 

!        89    formate  COMBUSTION  EFFICIENCY  (etac) •,fl2.3) 

j  write(*,*) 'Change  etac?  y-yes,  n-no • 

I  read (*, 20) yes 

I  if (yes.ge. 'y • .or.yes.eq. 'Y ' )  then 

j  write (*,*) 'Input  etac ' 

!  read(*,30)etac 

!  vol=aport*l  • 

j  go  to  99        . . 

j  endif 

!       c  calculate  the  nozzle  exit  area 

'  ae=ath*oex(2) 

c  calculate  the  required  perimeter  factor 

write (*, 19 )pc 
19  format (Ix, "Input  the  C*  with  Pcnom=' ,flO. 1, '  (F10.4) •) 

read(*, *) cstar 

xx=(pc*ath*3  2.17)/(cstar*rhop*l*3.14159*dpe) 
18  write(*,*) 'Input  the  estimated  required  web  (INCHES)  (F10.4) 

write (*,*) 'suggest  initial  guess  D/4 • 

(  read(*, *) webrqd 

r=webrqd/tb 

pfac=xx/r 

write (*, 21) r 

j        21    format (Ix, 'The  required  burn  rate  =',f6.4) 
!  write (*, 22 )pfac 

22  format (Ix, 'The  required  perimeter  factor=' , f6.3) 
write (*, 25) XX 

25    format(lx, 'pfac=',f5.4, '/r') 

write (*, 24 )tb  . 

24    format ( ix, 'tb  =',f7.2) 

write (*,*) 'Are  these  acceptable?  y-yes  n-no • 

read (*, 20) yes 

if (yes.ge. 'y' .or.yes.eq. 'Y' )  then 
pwet=pfac*pi*dpe 
wf=webrqd/(d/2) 
j=ath/aport 
ld=l/d 

go  to  23  i, 

else  '  ^     " 

go  to  18  ?    * 

endif  -  *   - 

20    format(al)  '■  *  .    . 

30    format(fl2.6)  ' 

23  write(*,40)cfc 
40     format (ix, 'THRUST  COEFFICIENT •,F12.3) 

write (*, 50 )etacf 
50    format(lx,  'ETACF.". •,F12.3) 

write (*, 60) isp 
60    format (ix, 'SPECIFIC  IMPULSE,   LBF  SEC/LBM •,fl2.3) 

write (*, 70) m 
70'    format (ix, 'MASS  FLOW  RATE,     LBM/SEC ',fl2.3) 

write(*,80)vp 
80    format  (Ix, 'PROPELLANT  VOLUME,   CUIN .',fl2.1) 

write (*, 10 0)ath 
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100 

ao5 

110 

120 

130 

140 

150 

90 

190 

170 

180 

200 


format ( 
write (* 
format ( 
write (* 
format ( 
write (* 
formate 
write (* 
format ( 
write (* 
format ( 
write (* 
format ( 
write (* 
format ( 
write (* 
format ( 
write (* 
format ( 
write (* 
format ( 
write (* 
formate 
return 
end 


IX, 'NOZZLE  THROAT  AREA,  SQ  IN •,fl2.3) 

,105)ae 

Ix, 'NOZZLE  EXIT  AREA,    SQ  IN ',fl2.3) 

,110)dpe 

IX, 'EQUIVALENT  PORT  DIAMETER,  IN ',£12. 3) 

,120)aport 

Ix,  'PORT  AREA,  SQ  IN  ',fl2.3) 

,130)pfac 

Ix,  '  PERIMETER  FACTOR ' ,  f  12 .3) 

, 140)webrqd 

Ix,  'REQUIRED  WEB,  IN -',fl2.3) 

,150)wf 

ix,  'REQUIRED  WEB  FRACTION ',fl2.3) 

,90)vl 

IX,  'VOLUMETRIC  LOADING ',fl2.3) 

,190)ld 

Ix,  'LENGTH   TO   DIAl-IETER  RATIO ',fl2.3) 

,170)tres 

Ix,  'RESIDENCE  TIME,  MSEC ',fl2.3) 

,180)j 

IX, 'THROAT  TO  PORT  AREA  RATIO  (J  FACTOR) . . . ' ,fl2.3) 

,200)r 

ix,  'REQUIRED  BURNING  RATE,    IN/SEC ',fl2.3) 


i 
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subroutine  loscf (etacf,tb, ath,pc, alpha, thetaex,erosr, 
*mf,nozle) 
c* ********************************** ************************ 

c  * 

C-  THIS  SUBROUTINE  CALCULATE  THE  NOZZLE  LOSS  COEFFICIENT     *  . 

C  (THRUST  COEFFICIENT)  USING  THE  METHODS  PROPOSED  BY  AGARD  * 

C  PROPULSION  AND  ENERGETICS  PANEL  FOR  METALLIZED  PRO-       * 

c  PELLANTS  (Reference:  Performance  of  Rocket  with  Metal-  * 

c  lized  propellants,  AGARD-AR-230)  * 

c  *         ■ 

c****** ******************************************* ********** 

c  Defination  of  varibles  in  this  subroutine 

nozzle  loss  coefficient 

percent  divergence  loss  » 

percent  kinetics  loss 

percent  boundary  layer  loss 

percent  two-phase  loss 

persent  submergence  loss 

percent  nozzle  erosion  loss 

nozzle  entrance/  nozzle  throat  area 

submergence  length/  length  of  internal  motor 

mole  fraction  of  condensed  phase  in  mole/lOOgrm  - 
i,c4,c6  : 

constant  in  emperical  relation 

diameter  of  particle 

chamber  pressure 

Isp  of  frozen  equilibrium 

Isp  of  shifting  equilibrium 

steel  nozzle  =  1,  ordinary  nozzle  =0 

burning  time  of  motor 

throat  diameter 

nozzle  expansion  ratio 

Isp  with  ideal  (at  design  point)  nozzle  expansion  ratio 

Isp  with  mean  nozzle  expansion  ratio 

integer  expnu,expnl 

real  etacf ,div,kin,bl, tp,sub,eros,astar,s,mf ,dp,pc, isps, ispm, ispf , 
+dt , ispu , ispl , erosr 

common  /scratc/  plot (5, 100) 

common  /socf/  spi(2) ,cf (2) ,oex(2) 
c  calculate  the  %  of  nozzle  divergence  loss 

div  =  50  *  (  1  -  cos  ((alpha  +  thetaex)  /  2  )) 
c  ispf  is  the  frozen  and  isps  is  the  shifting  specific  impulse  from        ' 
c  micropep  run  'at  nominal  pc 
c  Transfer  specific  impulse  data  from  MICROPEP." 

ispf=spi(l)  '^ 

isps=spi(2)  .  k 

c  Calculate  -  the  kinetic  loss 

kin  -  33.3  *  (1  -  ispf/isps) 
c  calculate  %  boundary  layer  loss 
c  Steel  nozzle  (1) ,  Ordinary  nozzle  (0) .  Enter  0  or 

if  (nozle.eq.l)  then 


c 

etacf 

c 

div 

c 

kin 

c 

bl 

c 

tp 

c 

sub 

c 

eros 

c 

astar 

c 

s 

c 

mf 

c 

cl,c2 

c 

c 

dp 

c 

pc 

c 

ispf 

c 

isps 

c 

nozle 

c 

tb 

c 

dt 

c 

expn 

c 

ispi 

c 

ispm 

c 

c 

cl 

= 

0. 

00506 

c2 

= 

0. 

0 

else 

cl 

m 

0. 

00365 

c2 

s 

0. 

000937 

endif 
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c    expansion  ratio  of  the  nozzle  for  shifting  equilibrium  from 
c    micropep  with  exit  pressure=ambient  pressure  at  design  altitude 
c     read  (*,70)expn 
expn=oex(2) 
pi=3. 14159 
dt=sqrt (ath/pi) *2 

bl  «{cl*pc**0.8/(dt*0.2) )*(l+2*exp(-c2*pc**0.8*tb/(dt*0.2) ))*(l+0. 
$016*(expn-9.0) ) 
c  calculate  the  two-phase  flow  loss 

c  the  mole  fraction  of-  condensed  phase  is  come  from  "MICROPEP" 
c  mean  particle  diameter  calculation  by  using  empirical  relation 
dp  -  3.39*dt**0.4692 
if  (mf.lt. 0.09)  then 

mf  =  0.09  ' 

endif 

if (dt.lt. 1.0)  then 
c3=9.0 
C5=1.0 
c6=1.0 
else 

if(dt.le.2)  then 
c3=9.0 
c5=1.0 

c6=0.8  *    - 

i  •      else 

if(dp.lt.4)  then  .    . 

c3=13.4 
;  c5=0.B 

I         -  C6=0.8 

'  else 

if(dp.le.8)  then 
I  c3=10.2  - 

I  ■  C5=0.8 

C6=0.4 
else 

c3=7.58 
c5=0.8 
C6=0.33 
endif  - 

endif  .  '  "•    ■ 

endif  •  ■ 

endif  '    ■  . 

c4=0.5 

tp=c3*mf**c4*dp**c5/(pc**0.15*expn**0.08*dt**c6) 
c  calculate  the  submergence  loss (NOT  USED  IN  PRELIMINARY  WORK  !1) 
c  the  nozzle  is  assumed  to  be  of  external  design  with  no  submergence 
c  loss.  s=length  of  submergence  /  length  *^  of  internal  motor, 
c     astar=(pi*d**2/4)/ath 

c     sub=0.0684*(pc*expn/astar) **0.8*s**0.4/dt**0.2 
C     default  to  0.0 

sub=0.0 
c  calculate  the  nozzle  erosion  loss 

c  calculate  the  ispm  (specific  impulse  at  mean  expansion  ratio) 
df=dt+erosr*tb*2 
dbar=(dt+df)/2 
expnm=(dt**2/dbar**2) *expn 
expnl=int ( expnm) 
c  Transfer  expansion  ratio  of  nozzle  and  specific  impulse  from  MICROPEP. 
ispl=plot ( 4 , expnl ) 
expnu=int ( expnm) + 1 
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ispu=plot (4 , expnu) 

ispm=(expnin-expnl)/(expnu-expnl)  *  (ispu-ispl)+ispl 
erosion  loss  calculation 

eros=(l-ispm/isps) *100 

write (*, 100) div 
100  format ( Ix, • PERCENT  DIVERGENCE  LOSS 

write (*, 110) tp 
110  format (ix, • PERCENT  TWO  PHASE  FLOW  LOSS 

write(*,120)bl 
120  format (ix, 'PERCENT  BOUNDARY  LAYER  LOSS 

write (*, 130) kin 
130  format  (Ix,  'PERCENT  KINETIC  LOSS , 

write ( * , 14  0 ) sub 
140  format(lx, 'PERCENT  SUBMERGENCE  LOSS 

write (*, 150) eros 
150  format ( Ix , • PERCENT  EROSION  LOSS , 

etacf=l-(div+tp+bl+kin+sub+eros)/100.0 

write (*, 160) etacf 
160  format  (Ix, 'THRUST  COEFFICIENT  EFFICIENCY 

return 

end 


',f6 

.3) 

•,f6 

.3) 

•,f6 

.3) 

•,f6 

.3) 

•,f6 

3) 

3) 

a 
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subroutine   input2 (pf ac , aport , web , rr5 , Igrain , tb , xx) 

c.  * 

C  THE  PRIMARY  PURPOSE  OF  THIS  SUBROUTINE  IS  TO  ES-  * 
C  TABLISH  THE  DATA  FILE  FOR  THE  CSD  "GRAINS"  PROGRAM  * 
c  * 

c* *********************************** *********************** 

y  character*80  tapeS 

!  Character* 1  yes 

real  Igrain, lspoke,kspoke,RR5 
integer  nmax , nspoke , ibb , iprnt 

write (*,*) 'The  following  will  be  the  data  input  for  • 
write (*,*) 'spoked  star  grain  designs  with  radial  sides  • 
open ( 4 , f ile= ' grainin . dat ' , status= ' unknown • )  • 

write(*,*) 'Run  DESIGN-1  mode  or  BURNBACK-0?  Type  1  or  0...  • 
read(*,*)i 
if (i.eq.l)then 
11      write (*, 40) web 
14      r=web/tb 
pfac=xx/r 
write (*, 13 )r 

13      format (Ix, 'Required  burning  rate,  IN/SEC •,fl2.3) 

write(*,50)pfac 

write(*,  *) 'Do  you  want  to  change  the  web?...y/n • 

read(*,5)yes 

if (yes.eq. 'Y' .or.yes.eq. 'y • ) then 

write (*,*) 'Input  the  required  web,  IN • 

read ( * , * ) web 
go  to  14 
endif 

write(*,*) 'Max  number  of  spokes?.  (13) ' 

read ( * , * ) nmax 

write(*,*) 'Initial  n\imber  of  spokes?.  (13) • 

read (*,*) nspoke 

ibb=l 

iprnt=0 

write (4, 10) nmax, nspoke, ibb, iprnt 

write (*, 30) aport 

30   format  (Ix, 'Required  port  area,  SQ  IN ' ,  f  12 

*.3) 

write (*, 40) web 

40   format  (Ix,  '  Initial  web  thickness  in  IN ' ,  f  12 

*.3) 

write (*, 50) pfac 

50   format(lx, 'Required  perimeter  factor •,fl2 

*.3)  , 

write (*,*) 'Input  Kspoke  -  Abstract  multiplier  on  ' 
write(*,*) •PHI2  -  if  KSPOKE=0.0  PHI2=0.0  where:' 

write(*,*) 'PHI2=ATAN(KSPOKE*WEB/RR3/COS(PHI3))    (0.0-1.0) ' 

read (*,*) kspoke 
write (*, 60) rr5 

60   format  (Ix, 'Grain  outer  radius,  INCHES •,fl2 

*.3) 

wr ite ( 4 , 2 0 ) aport , web , pfac , kspoke , rr5 

write (*,*) 'Input  initial  grain  inner  radius,  INCHES  (F7.3)  ' 

read(*,*)ril 

write (*,*) 'Input  initial  value  of  SRI  -  spoke' 

write(*,*) 'top  corner  radius  in  IN  (Default=0.25) • 

read(*,*)sril 

write (*,*) 'Input  initial  value  of  SR2  -  spoke' 
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write (*,*) 'base  coner  radius   in  IN    (Default=0.25) • 

read(*,*)sri2 

write(4,20)ril,sril,sri2 

write (*,*) 'Input  increment  on  SRI  in  IN  -  cannot  be  0.0...  • 

read(*,*)dell 

write (*,*) 'Input  increment  on  RRl  in  IN  -  negative  value..  ' 

read(*,*)del2 

write(*,*) 'Input  increment  on  NSPOKE  Default  "1.0 ' 

read(*,*)del3 

write (*,*) 'Input  increment  on  SR2  in  IN  -  cannot  be  0.0...  • 

read(*,*)del4 

write ( 4 , 2 0 ) dell , del2 , del3 , del4 

write (*,*) 'Input  tolerance  on  grain  angle  sum  total  -  %• 

write(*,*) 'of  360.(0.0-1.0) ..i • 

read(*, *)cirtol 

write (*,*) 'Input  tolerance  on  port  area  -  %  of  Aport  • 

write(*,*)  '  (0.0-1.0) ' 

read(*,*)pttol 

write (*,*) 'Input  tolerance  on  perimeter  -  %  of  perimeter..  ' 

write (*,*)  '  (0.0-1.0) • 

read(*, *)pertdl 

write ( 4 , 2 0 ) cirtol , pttol , pertol 

write (*,*) 'Input  Max  of  SRI  in  IN  .(F10.5) •' 

read(*,*)finl 

write(*,*)  'Input  Min  of  RRl  in  IN  .  (F10.5) » 

read(*, *) f in2 

write(*,*) 'Input  Max  of  SR2  in  IN  .  (F10.5) ' 

read(*,*)fin4 

write (*,*) 'Input  Max  of  web  in  IN  .(F10.5) • 

read(*,*) finweb 

write (*,*) 'Input  increment  for  web  in  IN.(F10.5) ' 

read ( * , * ) webstp 

write (4 , 20 ) finl , fin2 , fin4 , finweb , webstp 
else 

nmax=»0 

write(*,*) 'Number  of  spokes?. ..  (13) ' 

read(*,*)nspoke 

ibb=-l 

write(*,*) 'Do  you  want  optional  output?  y-yes,  n-no ' 

read(*,5)yes 

if (yes.eq. 'Y' .or.yes.eq. 'y') then 
iprnt=-l 

else  .,  . 

iprnt=»0  -^^ 

endif 

write(4, 10) nmax,nspoke, ibb,iprnt 

write(*,70)                  * 
70   format (Ix, 'Input  initial  port  area,  SQ  IN..(F10.5) ') 

read (*,*) aport  ■ 

write (*, 80) 
80   format (Ix, 'Input  initial  web  thickness,  IN.(F10.5) ') 

read (*,*) web 

write(*,90) 
90   format (Ix, 'Input  the  perimeter  factor. .  (FIO. 5) ') 

read(*, *)pfac 

kspoke=0.0 

write (*, 60) rr5 

write (4,20) aport , web , pf ac , kspoke , rr5 

write(*,*)  'Input  initial  RRl,  INCHES.  (FIO. 5) • 

read(*,*)ril 
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write (*,*)• Input  initial  value  of  SRI    (Default=0.25   IN.)...- 

read(*,*)sril 

write (*,*) 'Input  initial  value  of  SR2  (Default=0.25  IN.) 

read(*,*) sri2 

write ( 4 , 2 0 ) ril , sril , sri2 

dell=0.0 

del2»0.0 

del3=0.0 

del4=0.0 

write ( 4 , 2 0 ) dell , del2 , del3 , del4 

cirtol=0.0 

pttol=0.0 

pertol=0.0  I 

write (4,20) cirtol , pttol , pertol 

finl=0.0 

fin2=0.0 

fin4=0.0 

finweb=0.0 

webstp=0.0 

write (4, 20) finl,fin2,fin4,finweb,webstp 

write (*,*) 'Input  grain  length,  INCHES  (F10.5) • 

read(*, *) Igrain 

write(*,*)  'Input  initial  PHI2,  DEGREES  (F10.5) • 

read(*,*)dphi2 

write(*,*) 'Input  initial  PHI4,  DEGREES  (F10.5) • 

read  (*,*)dphi4 

write (4, 20) Igrain, dphi2,dphi4 

write (*,*) 'Input  #  of  steps  in  decimal  format' 

write(*,*) 'Input  #  of  burnback  steps  for  SR1>0 -• 

read(*,*)stepl 

write (*,*) 'Input  #  of  burnback  steps  for  PHI2>0  &  SR1=0....  ' 

read(*, *) step2 

write(*,*) 'Input  #  of  burnback  steps  for  LSPOKE=0  &  PHI2>0.  • 

read(*,*)st^p3 

write (*,*) 'Input  #  of  burnback  steps  for  • 

write(*,*) 'triangular  spoke  burnback • 

read(*,*)step4 

write (*,*) 'Input  #  of  burnback  steps  for  ' 

write  ( * ,  * )  •  remaining  fuel  burnout ' 

read(*,*) steps 

write  (4  ,20)  stepl,step2,step3,.step4,  steps 
endif  .,..  -    . 

5    format (al) 
10    format  (4i3)  •      '^        •     . 

20    format (5f 10. S)  >    ' 

close ( 4, status=' keep ' )  ii 

return     .  » 

end 
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